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ABSTRACT. ‘The approximate linear form of Einstein’s gravitational field is quantized 
by using an indefinite metric. It is shown that only two types of gravitons can be observed, 
though many more can exist in virtual states in the presence of interaction. 'The 
observable gravitons are shown to be particles of spin 2. Using the interaction 
representation, the interaction of the gravitational field with the matter field is briefly 
discussed. 


§1. INTRODUCTION 


HE present investigation is undertaken to carry out the quantization of 
| Einstein’s gravitational field, and to investigate the interaction of 
gravitational quanta or ‘gravitons’ with other elementary particles. For 
simplicity we shall first quantize the approximate linear form of Einstein’s 
field, while the exact treatment of the non-linear gravitational field will be given 
in a subsequent paper. Some work on the quantization of the gravitational 
field has been carried out in earlier days by Rosenfeld (1930), but since then 
very little progress has been made in this direction. 

For the quantization of the approximate linear gravitational field we shall 
follow the same treatment as has been applied to the radiation field in an 
earlier paper (Gupta 1950, to be referred to as A). ‘Thus we shall use an 
indefinite metric for the components of the gravitational field with a negative 
commutator, and then the gravitons corresponding to these components will 
be made unobservable by means of supplementary conditions. In this way 
we shall find that only two types of gravitons can be observed, though many 
more can exist in virtual states. Since the gravitational field has a large number 
of components, the present treatment is necessarily more involved than the 
treatment for the radiation field, but we shall not come across any new difficulty. 

It must be observed that the supplementary conditions play a vital role in 
our treatment. For, if there were no supplementary conditions, the present 
theory would involve observable states with negative probabilities, which would 
be physically meaningless. It appears, therefore, that the role of the 
supplementary conditions (or the coordinate conditions) is more fundamental 
than was originally intended by Einstein. 
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§2. LINEAR APPROXIMATION FOR THE GRAVITATIONAL FIELD 
The fundamental equation for Einstein’s gravitational field is given by 
Ri ho = — Ger ee (1) 


where the symbols g“’, R“’, and R have the usual meaning, « is a constant, 

and 7”” is the energy momentum tensor of the ‘matter’ field, where matter 

includes everything except the gravitational field. Since the covariant 

divergence of the left-hand side of (1) is known to vanish, it follows that the 
covariant divergence of T“” will-also vanish, whence one obtains 

‘ oe; 

Ox” 


where U* = T**(—g)? is the energy momentum tensor density for the matter 
field. Also, in order to express the conservation of energy and momentum, we 


write Gy to tae 
5,0 (Ent th) =05 Oe e- (3) 


ze qe Bes =0, 


ox" 


here t"” is the energy momentum pseudo-tensor density for the gravitational 
field, which satisfies the equation 
ot, = 1qrae Ou 
Fs —3@ 5 I asl ek sigh 8 s (4) 
Following Einstein (1918), we now obtain the linear approximation for the 
above field equations by putting 


Sasi ta 9. a ee (5) 


where «,, are the constant Minkowskian values for g,,, and the second and 
higher powers of « are to be neglected. Using the usual flat space-time tensor 
notation*, we can then write eqns. (1) and (4) as 


2 OTs Phy, Ohya, O)io 2 
Cees Ga a an,ox, +H" eats ~O hia) KT yy) cress (6) 
CLEP _ KOh,, $ 
Meare mc Sop (7) 


According to Einstein, we can also choose coordinate conditions or supplementary 
conditions given by 
Oh,» 1 Ohj; 


an, boxe eaiQyyicg ie OP Mest, Gite. Sea (8) 
so that (6) reduces to 
Gl i 38» [ Phi, =«T,,,. cee e ee (9) 
It is Convenient to put 
Piny =Yuv— 3 9 uv Vas eee (10) 
and thus we can write (9), (8) and (7) as 
i i 2 Sr i Secs (11) 
Ov 
oe an (tea: OT “egy vn idee (12) 
6). vo K OY y3 K Oya; " 
dx, 20,” Cae thy all ta en race's (13) 


* We take the space-time coordinates as x1, %2, V3, it. 
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It follows from the above equations that 
Lf yipPvig 10yrAVo0 1 5 (Fs Ovig 1 Oya, et) |. pee (14) 


t= eee 
ee ere ON, ON a OP Oy eae 2) Ox. ae 
In particular, we obtain for the Hamiltonian density 


qui OV 29 PVig 1 Oy 12 OY 0 Gee ee) | tee: (15) 


DB leGin cor 2-Or of 2\ 0%, One 2.0%, Oxe 


The eqns. (11), (12) and (13), describing a linear gravitational field, are so 
simple that it may seem surprising that they have to be regarded only as an 
approximation to a more involved theory. It should, therefore, be noted that 
the supplementary conditions (12) are compatible with the field eqn. (11) only 
in an approximate sense, for (11) and (12) give respectively 

2O%uv, OL wy 

Ne, Ox, 

which agree only in the first approximation. On the other hand, we shall show 

in the subsequent paper that for the exact non-linear gravitational field the 
supplementary conditions are exactly compatible with the field equations. 


OV» 
2 va gs 
andre! Ba 0, 


§3. QUANTIZATION OF THE GRAVITATIONAL FIELD 
Lagrangian density for the linear gravitational field may be taken as 
La 5 (Fede — 55h), Lier (16) 
4\ 0x, Ox, 20x,0x, 
where y,, and y are to be treated as independent variables. ‘This gives in the: 
usual way for the field equations and the Hamiltonian density 


Bat Fa 0, ary = 0, yee comers CE) 
a} [PB BY 1 (Pre Ptr 1 By By 
anc Say as AR Oi ON Ox, Ox, DOR ax ele ae ue) 


In order that (18) may agree with (15), we shall ensure by means of a 
supplementary condition, to be discussed in the next section, that the expectation 
value of y is equal to that of y,,... But for the purpose of quantization we regard y 
as an independent field quantity, so that we may easily split the contributions 
of the various components in the Hamiltonian (18). 

Now we have to obtain the commutation relations for y,, and y, keeping 
in view the fact that y,,=y,,. For this we note that the canonical conjugate 
of 4; is given by ere sl tvia * 

Ggtee aia ioe 2 Ge com) 
Therefore the commutation relation [y,,(r,t), a,(r’,t)]=2(r—r’) gives 
(yar, t), (A71,/0t) (r’, t)] =225(r —r’); or, using Schwinger’s notation (Schwinger 
1948), we have 


[yin(%), vir(%’)]=2eD(w— a"), ne ee (20) 
On the other hand, since y1, =21, we get 
Oot yA og “ee (21) 


ESC GAIED i) ae 
so that the commutation relation in this case is given by 


[yse(%), Yro(*’)]=tD(x— 4’), ne ww we (22) 
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According to (20) and (22), the commutation relations for y,, may be written as 
[Yul)s Yal%')] =4(8 18 90+ 809 va)D(H—%'). ve rae (23) 
Further, the canonical conjugate of y is given by 
2 leiilsaaal WG 5 
~ Weyid ~~ 4 Be 
and hence we obtain 
by(x), ye) |= = De) ae (24) 
Following A, we now regard y,, and y as self-adjoint quantities, and expand 


them as : 
Vuy = (2k)-U?{a,,,(k) eX" + al (k) ea eee a ae (25) 
k 


y= D2(2K)- UP fa(k) ek atk) eM (26) 
k 


where a dagger denotes an adjoint. Substituting the above Fourier expansions 
in (23) and (24) respectively, we get 
[ay (Kk), aho(k)] =8,184. + 8,105 va attoasntade 
fatk);al(k)i=—ibs-) 2m SRO) Vigieeer (28) 
Again, substituting (25) and (26) in (18), using the commutation relations (27) 


and (28), and omitting the zero-point energy, we obtain for the Hamiltonian of 
the gravitational field 


I Hav =Sk{4at,(k)a,,(k)—al(k)a(k)}. eee (29) 

k 

For the present purpose it is convenient to replace the operators a,,(k), 

Az9(k), dg3(k) and ago(k) by a7,(k), az2(k), a33(k) and ago(k), which are given by 

i(k) = ${441(Kk) — ao9(k)}, a3(k) = ue i(k) + ayo(k)}, 

1 

A33(K) = W/2 A33(k), Ayo(K) = 5 Aqo(k). 

“according to (27) and (30), the commutation relations for the new operators are 


Caiu(k), axl (K)=1, — Laa(k), ai(k)]=1, | 
Lajs(k), aft (K)]=1, Caroll), affk=1, fo oe 


and the Hamiltonian (29) may be written as 
| Hedo =EKfa (ke) aya( ke) + ak) ag(e) + 4k) x(k) 

—alo(k )ayo(k) — a3o(k)aoo(Kk) — ao(k)agq(k) 

+ ayi(k)ay)(k z 433(K)azo(k) + a33(k)a53(k) 

+ apg(k)aoo(k)—al(k)a(k)}, nee (32) 
where we have used the fact that ae a,,(k). Thus we have, in all, eleven 
types of gravitons corresponding to the eleven independent components of 
Yu» and y. We shall refer to the gravitons corresponding to the operator a,,(k) 
as the a,.-gravitons, and so on. 

According to (27) and (28), the commutators involving aj(k) and a(k) have 
a negative sign. ‘Therefore, as explained in A, we have to use an indefinite 


metric for the components yj) and y of the gravitational field, while the remaining 
components are to be treated in the usual way. ‘Thus, all the operators a,,(k), 
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A a(K), a3x(K), ay0(Kk), Aao(K), ago(Ke), aii(K), azo(K), a33(K), apo(k) and a(k) appear 
as absorption operators, while their adjoints are all emission operators. It is 
to be noted that though now all the gravitons have positive energies, some states 
will occur with negative probabilities. For, due to the use of an indefinite 
metric, we have to normalize any state as 


Es 1 ae caet eae 20 REY JER (33) 


where 19, %z9, %39 and m are the numbers of the ajo-, dyg-, Ago- and a-gravitons 
respectively in the state Y’. This occurrence of negative probabilities, however, 
does not raise any difficulty of physical interpretation, for we shall show below 
that for all real states we have 14) =Mg9 =Ng) =n =0. 


§4. THE SUPPLEMENTARY CONDITIONS 


In the absence of interaction the supplementary conditions for the 
gravitational field may be taken as 
(7) 


OY ws 9, 


AAR (34) and (Vg =) iO) festa) 
0x, 


where y\;) and y+ are the positive frequency parts of y,, and y respectively. 
It follows from (34), (35) and their adjoint equations that <dy,,/0x,)=0, 
(Yun) =<y)- ‘Thus the supplementary conditions (34) agree with (12) in the 
classical limit, while (35) ensures that the Hamiltonian (18) is equivalent to (15). 

Substituting the Fourier expansions (25) in (34), and choosing the x3 axis 
along k, we obtain 


{d@g,(k) —ag,(k) PP =0, © so..02 (36a) 

{gal K) c= Gyg( KK) Y= 0, ca (36 5) 
{agq(K) — dgg(k) j= (y/2a5o(k) — doa(k)}F=0, — ...... (36 ¢) 
(430(K) — ago(k) j P= {ago(k) — V/2aoo(ke) P= 0... (36d) 


As shown in A, a solution of (36a) consists of a normalized state containing ne 
@z,- and a ,-gravitons, and a superposition of a series of redundant states with 
arbitrary coefficients and zero normalization. ‘These redundant states, however, 
do not contribute to any observable effect. Moreover, since the normalization 
is conserved in course of time, a redundant state will always remain redundant. 
Thus for all practical purposes we may ignore such states, and regard the 
@y,- and do,-gravitons as entirely absent in a non-interacting gravitational field. 
Similarly, (365) implies the absence of ag9- and d.-gravitons in the absence of 
interaction. The situation regarding eqns. (36c) and (36d) is slightly 
simpler, for both of these equations involve the operator dy,(k). It can be easily 
seen that the only states satisfying (36 c) and (36 d) are those which do not contain 
ANY Apg-, Aj3- and ajy-gravitons. 

We now consider the supplementary condition (35). Substituting the 
Fourier expansions (25) and (26) in (35), we get 


{a,,(k)—2a(k)}¥=0. = = sacee . (37) 
Using (30), (36c) and (36d), we obtain from (37) 
Ch Ds AG ae Lh de Spe rc (38) 


Again, as explained above in the case of eqn. (36a), (38) implies the absence 
of aj.- and a-gravitons in a pure gravitational field. 
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We have seen that the supplementary conditions (34) and (35) enable us to 
eliminate nine types of gravitons, so that for a pure gravitational field we may take 


if / / 
No a Ng — Noy — Nog eae Nog = Ngo — N32 = Noo =n= 0, sya aliouele (39) 


where 13 is the number of dy3-gravitons, and so on. Therefore, in a pure 
gravitational field, only a,,- and aj,-gravitons can exist. It follows that the 
state of vacuum for the gravitational field may be defined as that containing no 
@o- and aj,-gravitons, because the other types of gravitons will be excluded 
by the supplementary conditions. Moreover, since according to our treatment 
all the components of y‘!) and y+) contain the absorption operators, the state 
of vacuum will satisfy the relations 


YL =On 1 eee Medien 2 ee (40) 


$5. SPIN OF THE GRAVITONS 
In the last section it was shown that only two types of gravitons can exist 
in a pure gravitational field. We shall now discuss the spin of these gravitons. 
From the Lagrangian density (16), one can obtain the canonical energy 
momentum tensor #,, for the gravitational field as 


ke Oya OL dy OL 
lig= = ax, O(0y19/0%,) = ax, A(0y/dx,) FO ybss nate there (41) 
This energy momentum tensor, however, has to be modified according to the 
treatment of Belinfante (see Wentzel 1949). For this we consider an 


infinitesimal Lorentz transformation . 
dx,=S,,%, With da,,=—0W,,, sees (42) 
and find a tensor /,, , given by 
ae SY an ie OL 
vee OW, O(OY7/OXp) 


. Oi OV ux 
= Ox, Yon Ox, Vie ee Sg BFale antes nis (43) 


Then the modified energy momentum tensor 0,,, will be 


Le 

Oy =tyy— 5 Gn, Suro t Fon, ot Love). ae Me (44) 

and the angular momentum density is given by 

Mav, 0 = % 9 99 — XP ug 
= Xybyg— Xotuot tay, o- at a ealet on 
Thus we obtain for the components of the angular momentum 

P= —i| Mg,.do= Pt Ply eaceal 40) 
with raha a | do (x,T x4 —%nT a), He TAT 


, : : OVE, Oy; 
Pa= ~i| dv fa, 4=i | de (5 22 ~ yn 5H"), deca (48) 
4 4 
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where Pj, is the orbital angular momentum, and P%, is the contribution due 
to the spin. Substituting the Fourier expansions (25) in (48), we get 


Pip= DDS | dev [— {au(k) el"! 4 ah (k) eer} 
k k’ 


a 
x {a,,(k’) Cr es eat kK or Rete 
+ faz, (Ik) ee *-#9 4 at (ke) eter 
x {a;,(k’) on List Ade ai,(k’) eT tk 3 rey] 


= " : [a,,,(k)aj,,(k) iar ai,(k)a zn( K) 7 ax,(k)a},(k) a d.,(k)a,,(k)] 
ae ama Kye k) seal (Ria (k) 27 elmira ts Ale (49) 


Since the real gravitons correspond to the operators 
Q(k) and  aj,(k)=3{a11(k) — a9(k)}, 
we have to consider only that part of Pj, which involves the operators 
44(k), @,,(k), dy9(k) and their adjoints. ‘These operators occur only in the 
component P/,, and their contribution is given by 


Pig= i[a},(k)a,,(k) 7 Ayo( k)ay2(k) — a} ;(k)ay1(k) fos Ajo(k)ay9(k)] 
= 2ilaje(k)ay(k)—ayi(k)ay(k)], ne (50) 


In order to separate the contributions of the two independent components in 
the above equation, we introduce the operators 


as(k) = Sy fahi(k) tall}, a(k) = Sy fal) igh... (53) 


where, according to (27) and (31), we have 
[a,(k), a',(k)] =[4-(k), at (k)]=1, | sr (52) 
[a,(k), a‘ _(k)] =0. 
In terms of these new operators we obtain 
Pr» ai 2[al, (k)a,(k) — a! (k)a_(k)] 


ST ee ae ee ee (53) 


where 7 ,(k) and n_(k) give the numbers of the gravitons corresponding to the 


operators a,(k) and a_(k). 
This shows that the gravitons are particles of spin 2, and they have two 


independent spin states with the spin axis parallel or antiparallel to the direction 
of motion of the gravitons. 


§6. TREATMENT OF INTERACTION 
It is evident from (11) that the interaction of the gravitational field with 
the matter field may be taken into account by modifying the Lagrangian 
density (16) as 
K 
‘moe i(Z Pus _ 5 oe 5] SG REC Lae (54) 
which gives for the field equations 
CP yu =K Lys CPy =x, eae) 
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In the presence of interaction, the supplementary conditions (34) and (35) 
have also to be modified. Following A, we may state the supplementary 
conditions in the Heisenberg representation as 


qe | =o, Vag Vi Oe gerd (56) 
ra 


where [ ]+ denotes the positive frequency part. ‘These conditions are, 
of course, identical with (34) and (35) in the absence of interaction. But the 
necessity of writing them in the form (56) arises from the fact that in the 
Heisenberg representation y,, and y cannot be split into the positive and the 
negative frequency parts in the presence of interaction. However, such a 
splitting is still possible for dy,,/@x, and y,,—y. For, according to (55), 
Oy,,,/0x, and y,,,—y still satisfy the wave equation for plane waves: 


OY u» 
i ae ot [SI patent) =0, 


where we have taken into account the approximation used in §2. 

When we pass over from the Heisenberg representation to the interaction 
representation, the required supplementary conditions may be obtained in the 
same way as in the case of the radiation field. Following the treatment given 
for the radiation field (Gupta 1951)*, we easily find in the interaction 
representation 


(Me +k | x D“(x% — x')To,(x') ac’) Y(t) =0, 7 


Ox, Yt 
(Yun YP) P=, 


where D(x—«x’) is the positive frequency part of D(x—«x’). It should.be 
observed that, due to the occurrence of the 7y,-term in (57), those gravitons which 
could not exist in a pure gravitational field can now appear in virtual states in 
the presence of interaction. 

For practical purposes it is most convenient to use the invariant perturbation 
theory based on the interaction representation. For this we require, in addition 
to the commutation relations (23) and (24), the vacuum expectation values of 
the anticommutators. Using the vacuum conditions (40), we easily find, by 
following Schwinger (1949), 

({¥y%), Yaol&’)} Yo =F (5,29 vO tr Su va) D(x it x), | 

C{y(%), 1(2')} 0 = —4DO(x — x’). } 

It should be noted that the use of an indefinite metric does not make any 

difference to the calculation of the vacuum expectation values, for the 

indefinite metric operator 7, defined in A, is just equal to unity for the state of 
vacuum, which does not contain any gravitons. 


§7, LORENTZ-INVARIANCE 
In the present paper we have treated the components yo, of y,, by means of 
an indefinite metric, while the other components have been treated in the usual 
way. ‘I‘hus, though the expectation value of y,, is always real, we have to 
regard y,, as a complex tensor with yo, antiHermitian and the other components 
Hermitian. In this way the Lorentz-invariance of our treatment is not quite 


* The letter, referred to above, contains an error. There the relation preceding eqn. (9) should 
read [8D(~w—x’)/dt], y= — B(x —Xy’) B(x — XQ’) B(03 — 20’). 
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obvious. However, the Lorentz-invariance of the present theory can be easily 
established by splitting up y,, into the Hermitian and the antiHermitian 
components in a more general way, as we have done in A. 

For this we take y,, to be of the form 


Ypv=NONPAP + NON PAM 4+ NON GAG) — NON MAG 
+ NONMACD +4 NONOACY + NON@ACD— NON AGO 
+ NONPASD + NONOPAS + NONOAGB™— NON OAC 
—NONPAM— NONE — NONAG4 NONMAOM, (59) 
where, as in A, NY, N, N® and N are a set of orthogonal unit four-vectors, 


“2 

Ree Chon being space-like, and N° time-like. Further, A@?, AQ®, . . 
are scalars, of which A®®, A@%, AG), AOD, A) and A®) are antiHermitian, while 
the rest are Hermitian. Thus, according to (59), we have split y,,, into sixteen 
components in a formally Lorentz-invariant way. We can now follow the 
entire treatment of the present paper with two differences: Firstly, instead 
of y,,, we have to expand the A’s into Fourier expansions of the form (25). 
Secondly, instead of taking the x, axis along k, we have to choose, as in A, 
NY, N®, N© in such a way that 

NE GROIN ahead) Mile Note RRNA ser Dy tae (60) 
The treatment for y, of course, remains unchanged, because y is a scalar 
independent of y,,. In this way the vacuum conditions (40) and all the other 
results remain unchanged. 

Since (59) is a tensor equation, all our results will now be formally 
Lorentz-invariant. Still, in order to establish the relativistic invariance of our 
treatment one further point remains to be clarified. According to (59), if we 
choose a frame of reference in which N{° =(0, 0, 0, 7), ie. the time-axis 
coincides with N“”, the components yo; of y,, will be antiHermitian while 
the other components will be Hermitian. Since this property of y,,, is not 
preserved in all frames of reference, our treatment has singled out those frames 
of reference in which the time-axis coincides with N°. Thus, in the present 
theory N<° appears as a preferred direction. However, this does not matter, 
for N{° can be chosen in an arbitrary way, and the results of physical interest 
are independent of its direction. In fact, after the vacuum conditions (40) have 
been established, N° never occurs explicitly in the interaction representation. 
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ABSTRACT. A new coordinate condition has been considered and the equations of 
motion of two bodies in their weak and quasistatic gravitational field have been derived 
to the second approximation. A number of additional second order terms appear in 
the new equations which, however, have no influence on either the secular motion of the 
centre of mass or the rotation of periastron of the two-body system. 


$1. 
N recent years two independent methods have been developed for the 
| derivation of the equations of motion of bodies in weak and quasistatic 
gravitational fields from the field equations of general relativity. The 
method of Einstein and his co-workers (1938, 1940, 1949) considers the bodies 
to be spherically symmetric singularities, and derives their equations of motion 
from the field equations of empty space surrounding the singularities. Fock 
(1939) and Papapetrou (1951a) employ the field equations inside the bodies 
as well as in empty space; their method can, in principle, be applied also to 
bodies without spherical symmetry. 
An essential step in both methods is the use of four first order differential 
equations which are subsequently referred to as the coordinate condition. The 
coordinate condition used by Einstein and co-workers is 


hg yh) pat ee (1.1) 


where, as usual, h,,=g,,—7,, and an index after a comma denotes ordinary 
differentiation. Fock and Papapetrou use the coordinate condition 


gh) a0), anisole PL Seer ae (1.2) 


Both coordinate conditions are identical to the first approximation, i.e. when 
retaining terms up to the order 1/c? only; however, they differ in the higher 
approximations, i.e. when terms of the order 1/c*, 1/c® etc. are retained. 

The question arises whether the choice of the coordinate condition will 
influence the form of the equations of motion. Einstein and Infeld (1949) 
have proved that if the coordinate condition (1.1) is changed arbitrarily, but in 
the second approximation only (i.e. being kept unchanged in the first approxi- 
mation), ther: the equations of motion up to the second approximation will 
remain unchanged. ‘This result explains why the methods of Einstein and 
Fock lead to the same equations of motion, though the conditions (1.1) and 
(1.2) coincide in the first approximation only. 

A systematic discussion of the role of the coordinate condition has been given 
by Papapetrou (1951b). He finds that in the first approximation the equations 
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of motion do not depend on the coordinate condition, the assumption of a weak 
and quasistatic field leading necessarily to the Newtonian equations of motion. 
To any higher approximation the equations of motion will depend on the 
coordinate condition, though for the equations of motion up to the th 
approximation only the form of the coordinate condition up to the (n—1)th 
approximation will be of importance. 

In the present paper a new coordinate condition, which differs from (1.1) 
and (1.2) already in the first approximation, will be considered and the corres- 
ponding change of the equations of motion to the second approximation will 
be calculated explicitly. The new coordinate condition is g’”,=0, or to the 
first approximation nh 


The equations of motion will. then, acquire a number of additional second order 
terms. ‘These terms will be calculated by means of an appropriate transformation 
of coordinates. 


2, 

We denote by an asterisk a system of coordinates in which the coordinate 
condition (1.1) holds. This is the system employed by Einstein, Infeld and 
Hoffman (1938) in which the equations of motion of two bodies up to the second 
approximation have been calculated. With a slight change of notation we write 


Ajo =hiot hiot +++: 

n= Hon tlton "Hees ee NE Sa eee (2.1) 
Os peat RE * 

high PBS et. | 


mn on 


the subscripts under the terms on the right-hand side indicating the order of 
magnitude of each term.t 
We now consider a transformation of the form 


aera seas oe ang Be (2.2) 
r=1" 
which differs from the identical transformation x*“=x" by terms of at least the 
first order. It will be shown in §3 that, if we retain the basic assumption of the 
field being weak and quasistatic, a transformation of the type (2.2) will be 
sufficient for fulfilling, in the frame x“, the new coordinate condition (1.3). 
Moreover, we may retain Einstein’s demand that quantities with an odd (or even) 
number of indices 0 be developed in an odd (or even) power series. ‘Thus our 
transformation will be of the form 
xt = x94 G94 g9+...., | 
oes wer 23) 


ao*® = we + a® + ae + Saas j 


It is seen immediately that in order to get the equations of motion in the 
frame x", to the second approximation, it will be necessary to calculate the 
quantities g° and a’ only. This calculation will be given in §3. We simplify 


+ As usual, Latin indices take the values 1, 2, 3 only, while Greek indices take the values 1, 2, 3 
and 0 
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the present calculations by anticipating a result of §3, namely that a°=0. Thus 
aii ha Ret en Maeda AM ee ar (2.3 a) 


According to Einstein, Infeld and Hoffman (1938) the equations of motion 
a he 1 ; ee 
of a body of gravitational mass m and space coordinates 7** in the presence of a 
second body of mass m and coordinates ¢** are 
dy _g @ (1) _g [[antednts | 3dbtedte ,dytedi" ym _ ) 
den haat ep || ats an? Ueda eden dee ane Ae 


; (2) tush a a ann di ten de 


‘ oT dt* \ dt” dt* ae are ain 
Bich Aik eta fas seupal EGE 
an*s \r* 55 Sakae ®t ee 2.4 
x on*8 (=) 9 20n*”"dn*80n*" dt* dt* | ) ( ) 
where (GF =(n*8— £*9)(n** — C8), 


The equations of motion of the second body are obtained by replacing m, ™, nee 
and ¢*§ by m, m, *® and n*s respectively. 

The equations of motion in the frame x“ will be obtained by introducing 
(2.3) into (2.4). We notice that the left-hand side of (2.4) is of the order 1/c? 


since t* actually stands for ct* and m, m are the masses multiplied by G/c?), 
while the right-hand side is of the order 1/c*. Since eqn. (2.4) is correct only 
up to terms of the order 1/c* we see that the transformation (2.3) will change 
the left-hand side only. Let us write 
5y° ="* — 4° =value of g* for x8 =7°, { 
dfs = 6*8 — (8 = value of g° for. x= 095) 


and fi 22(7° — (i=) | 7 eee (2.6) 
An elementary calculation gives 


donee dn” da” F 
ae ~ ae + geen"): 


O/T 02 (AN 508 NOD a mee 
ar (r=) ~ ay 7) + agra (7) 4 -*)~ a7) age 


Hence in the frame x“ we shall have in the left-hand side of the equations of 
motion the additional terms 


d? Dale Oteny (al 2 0 /1\ a 
On = ee s__ S7s pata) seat) eet, 8 
7p | n”) ma noi @) (88 — 58) + ice (7) On (CS) Ree (2.7) 
We notice that there has been no necessity for assuming any change of the 
masses m and m in the transformation from x*" to «“. This is why we did not 
put asterisks on m and m in eqn. (2.4). 


+ The more general assumption would be that in the transformation from x*# to x“ described 
by (2.3) we should have m*=m-+ 8m where 8m is small of the first order. Hence a term 


5m 0/0”)(1/r) would appear on the right-hand side of the final equations of motion (3.14). Sucha 
term might have some physical meaning only if it would be possible to determine the quantity m 
by some independent method. Actually the mass of a star is determined by observations of the 
motion of other bodies in its gravitational field, and is by definition the coefficient of the term 
(0/én”)(1/r) in the equations of motion. Thus the splitting of this coefficient into m-+ 8m will 
remain arbitrary and, therefore, physically meaningless. 
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§ 3. 
The functions Q and q will be determined from the demand that in the 


frame «“ the new coordinate condition (1.3) be fulfilled. Using the transformation 
formula connecting g,,, with gi, we find, 


— B* 0 ; 

fron = hoo + 22° 6 
= 7, 0 

hon =hon+4 ee mee? a ae hee ee (3.1) 
S- Doe Ms ot 0 0 0 * 0 * 

Ann himn “ 7” ‘ mt DD mt mhon te nom: 


The values of the quantities h/, are previously known to be (Einstein, Infeld and 
Hoffman 1938), 


Po 2 
m * * 
hog = =e iP ig ? hy, =9, Bnn = Smnbtoo ee AO PG. (3.2) 
where Se maa ay (rae el (xe) (3.3) 
The coordinate condition (1.3), up to terms of the order 1/c?, is 
hos, s=9, Fino, 0 — Bns,s =9- cee eee (3.4) 


Introducing (3.1) and (3.2) into (3.4) we find that q° and g* must satisfy the 
differential equations 


Omer n Ss 0 0 0 sai hes 
Ay ss =9, we st Zen ta ,on—(2 a ,s),s = Hoo, »: FO O00 (3.5) 


The coordinate condition is valid throughout space, i.e. inside as well as outside 
the material bodies. Therefore it follows from the first of eqns. (3.5) that a =) 
The second of eqns. (3.5) then becomes 


QB" sot @en=hon- tte (3.6) 


a 


If eqn. (3.6) is differentiated with respect to x", we find 2a”, nss = Moo, ane sinee 
this equation is valid throughout space, it follows Go = doo: Equation (3.6) 


then reduces to : 
Q” 3s = Thoo, PO pel i eee ee (3.7) 
The solution of this equation can be written in the form 
dw 
air 6 a? a eer text (3.8) 
where the function w satisfies the equation 
®, os = doo Aaa (3.9) 


To solve this equation we must first of all express hoo in terms of the 


coordinates x". Remembering (2.3) we see that this is done by simply neglecting 
the asterisks on the right-hand side of the first equation in (3.2) and in eqns. (3.3). 
Now oo is the sum of two terms and (3.9) is a linear equation. Therefore, w 


will also be the sum of two terms, one for each of the two bodies. 
Let us, then, first discuss eqn. (3.9) for one body only. In order to determine 
w we shall have to consider the values of Ago throughout space (i.e. outside as well 


as inside the body). Temporarily, we denote by 7 the distance from the centre 
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of the body and by p(r) the mass density at the distance r. From the equation 
Ahoy =87p, we find the general formula 


phy | pidr—8n{prdr eases (3.10) 
a ’r Jo r : 
Since hoo is a function of r only, the same will be true also for w. Remembering 
et Pu 2de 
W 9s= Aw = LED + Fi eS) 
we find from (3.9) and (3.10) 
@w 2dw 


ae ies = [prt dr—8n | “pr dr. 


It will be sufficient to integrate this equation once only since 


i.e. for the calculation of g” only the quantities dw/dr are needed. 

The constant of integration will be determined by the demand that dw/dr be 
finite throughout space, since otherwise the transformation (2.3) would have no 
meaning. We find, 


dw _ Px als LTT ct, 4ar (° 
= —2n | pr drt 5 pr dr arden Gates (3.11) 
It follows from this formula firstly that 
dwjdr=O— tor: r=0l- 2 sae) Pe aoe (3.12) 


Outside the body the last term of (3.11) vanishes. For r>R, R denoting the 
radius of the body, the second term on the right-hand side of (3.11) becomes 
-R 


negligible compared withthe first. ‘Therefore, remembering that 47 | pr? dr=m, 
“0 


dujdr=— tm» for’ irSRy Sha eee (3.13) 


The total w will now be w=w+0, w being a function of r and & a function of 7. 
Remembering (3.8) and taking into account (3.12) and (3.13) we find for the 
quantities 5y* and 8f* defined by (2.5) the values 6n%=— Lin (ns — &)/r, 
865 = 4m (78 — ¢*)/r where r now has the value defined by (2.6). 

From these expressions we calculate the additional terms (2.7) which appear 
on the left-hand side of the equations of motion. In calculating these terms we 
substitute the Newtonian approximation for d?n"/dt? and d?¢"/dl*, i.e. 

OTe 


SLO i! Ge ote LNG Wee’ all 
at On" \r)? db =i apn ; ~ are ah 


Sr 


Taking into account the relations 


20) oe eee 
OnP XT Ons dn” — On ORENG On? = rn” =) 


we find that the additional terms (2.7) have the value 


ea mr (dn* _ al) (der _ at 
ee 2 dn” On” Ons \ dt a i i). 


r dt at 
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The final equations of motion in a coordinate system which satisfies the 
coordinate condition (1.3) are therefore 


dn” a /l d 8 dts sds om on 
ee (=) =f ores See ea woah 
dt® ” On* dt dt ‘2 dt dt dt dt r On" \r 


ac* dy” dy” aye 4 46" dts 
+[4 qF -¢) Paar poe ia ae (3) 
o'r (3 al", ldnidn? dl? dn ‘ \ 


Sn" Ontdn’ \dt dt | 2 at di dt dt 


94: 

Robertson (1938) has integrated the equations of motion (2.4). Robertson’s 
method may be repeated for the eqns. (3.14). The result of this calculation, 
which need not be described here in any detail, is that the additional terms in 
the new equations of motion do not influence the two quantities which 
Robertson discusses, i.e. the secular motion of the centre of gravity and the 
periastron advance of a double star system. If one assumes that these quantities 
are the only observables of the motion, then our result would mean that the 
coordinate condition (1.3) is fully equivalent, from the physical point of view, to 
(1.1) or (1.2). Such a statement, however, is not plausible at all since one would 
expect that an observer, on one of the two stars, for example, might be able to 
observe many more details of the motion; if this is so, then—according to the point 
of view expressed by Papapetrou (1951 b)—the conditions (1.1) and (1.3) need 
not be equivalent. 
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ABSTRACT. A partially polarized beam of particles is represented by a statistical 
operator in spin space. For nucleon—nucleon scattering a theorem stated by Wolfenstein in 
1949 connecting the two simplest experiments involving polarized nucleons is shown to 
depend on the invariance of the interaction potential under time-reflection. 


of polarized beams and the polarization of scattered nucleons are 
logically the next step in the unravelling of the nucleon-nucleon 
scattering matrix. Such experiments have been proposed by Wolfenstein 


IN patria bea scattering experiments which study the scattering 
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(1949 b) and by Lepore (1950), who have evaluated the polarization produced 
in a single scattering, and the azimuthal asymmetry in scattering of a polarized 
beam, for some tensor—force interactions. Here it is our purpose to discuss the 
concept of polarization of a beam and to study the structure of the nucleon— 
nucleon scattering matrix. 

Since a wave function % will always describe a particle completely polarized 
in some direction, a partially polarized beam cannot be a pure state, but is 
in fact a statistical mixture of pure states. Such a mixture may be specified by 
giving the set {*} of normalized states occurring, together with the relative 
frequency P(«) with which each % is represented in the beam. If a basic set 
of states {/,} is chosen, so that we may write 4*=2,C,%,, the expectation 
value of an observable O in the beam is 


O=%,P(a)(p", Of") = ZijpP(a)C* C70, 
=2,,U;,,0,= Trace (UO), 0 OCR eee (1) 


PAM) 
where U is the matrix (U,,) where U;,=,P(«)C,;*C,". From the properties 
of the probability P(«), U,=Trace (U)=1. This statistical matrix U now 
contains all information concerning the beam necessary for prediction of the 
results of experiments with the beam. Except for the case of a completely 
polarized beam (which is a pure state), many different probability distributions 
P(«) and sets of states {%} will lead to the same matrix U. 

For example, a linear beam of spin 4 particles may be represented by a 
statistical operator U in spin space. Such a 2x2 matrix may always be written 
in the form A+B.o. Since the normalization condition implies A=}, it is 
convenient to write U=}(1+p.ce). The pseudo-vector p thus specifies 
completely the polarization nature of the beam and may be referred to as the 
polarization of the beam. At once p=Trace (sU). 

For spin greater than $, U will be more complicated, e.g. for deuterons 
U=({1+),S,+ T,,(S,S;+ S,S,)}/3, S; being the spin operator for total spin 1, 
so that the specification of the state of the beam requires both the pseudo-vector p;, 
and the symmetric traceless tensor 7;;. In particular, for such a beam it is 
possible that p,=0, yet 7;;40, so that the ordering in spin-space is generally 
more complicated than a simple polarization. ‘These considerations may 
clearly be generalized to any total spin, and it is noteworthy that the analysis 
is at once applicable also to relativistic particles, for example electrons. 

In spin space the scattering of a particle may be described by an operator 
of(6,k;,k;), which transforms the initial to the final spin state, k,, ky being unit 
vectors in the initial and final directions. ‘The outgoing beam ky will be 
represented now by the statistical matrix 


Uf" = (G7 Uo i tate (G2 en) aoa ee (2) 
so that, for example, the polarization after the collision will be 
p’ = Trace (oU’)=Trace (ao% Uc#")/'Trace (e% Uc¥!). _...... (3) 
On the other hand, the differential cross section for the collision is 
do/dQ. = Trace (cf Uc") =4{ Trace (e%o/") + p. Trace (o%ac¥")}. ...... (4) 


It is of interest to compare these considerations with earlier work on the 
scattering of polarized thermal neutrons by unpolarized nuclei (e.g. Borowitz 
and Hamermesh 1948) or by ferromagnets (e.g. Halpern and Johnston 1939). 
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In the former case the scattering operator is o&% =a+ba.l, | being the nucleon 
spin operator. By evaluation of (3) for unpolarized nuclei it is easily shown 
that the direction of the neutron polarization is unchanged by the collision, 
though the magnitude of the polarization is always reduced. A choice of basic 
spin states directed along this constant direction leads to the representation 
usual in this study, i.e. by a polarization factor (P,,— P_)/(P,.+P_), P,, P_ being 
the probabilities of the two spin orientations. A similar representation has been 
used in the study of scattering by a ferromagnet, in which case the justification 
lies in the precession of the spin about the internal magnetic field by a random 
amount between collisions, as noted by Halpern and Johnston (1939). 

The elastic collision of two nucleons is of special interest. The operator 
describing the scattering from relative momentum direction k, to direction k, 
may be written e/(6,, 6,; k;, k;). Since the Hamiltonian will be invariant for 
rotations and reflections of the coordinate system, the most general formf possible 
for of is 


oF =A + Bo, .k; x k-+ Co,.k; x k,-+ Do, . ko, .k; + Eo, . ko, . ky 
a Fo,.k,o,.k-+ Go, .k,o,.k;+ Ho, .k; x k,o,.k; x kp+Jo,.65, atodelenate (5) 


where A, B, C, D, E, F, G, H, J are functions of k,.k;=cos 6. 

It is natural to suppose that the Hamiltonian is invariant also for reversal 
of the time axis, which implies further restriction on the structure of e%. Under 
time reflection, e” behaves as follows: 


t+t' = —t : (04, 65; k;, kp) > 0%’ =e/(—0,, —o,; —k;, —ke). ...... (6) 


However, the time-reversed scattering process is the scattering of an incident 
wave of direction —k, to final direction —k,, so that 


P(— 01, — G2; — ki, — ks) =07(6,,6,; —ke, —Kki). ss... (7) 
This equation implies that D= EF, F=G in the form (5), so that 
of = {A U(o, + 6,)-k; x k-+ V(o, = o,)-k; x kp+ D(o,.k;o,. k;+ Sy e ko. k;) 
+ F(o,.kjo,.Kp+0,-kpo,.k;)+Jo,-0,}, = = 9 «2.20, (8) 
where U=(B+C)/2, V=(B-C)/2. 

It has been stated, but not proved, by Wolfenstein (1949a,b) that the 
polarization produced by scattering of an unpolarized nucleon is equal to the 
proportional value of the term not azimuthally symmetric in the differential 
cross section for scattering of a completely polarized beam by unpolarized 
target nucleons. This important connection between the two simplest 
polarization experiments is, from (3) and (4), equivalent to the statement that 
Trace (6,0%0/") = Trace (o%a,0/"). By direct calculation, using (8), each of 
these quantities is 
Trace (o%o,0%") = Trace (6,07 0%") =k; x kp 82{A"(U+ V) + J"(U—V)}...(9) 
It is significant to note, however, that this relation is not valid for the more 
general scattering operator (5). This connection between the simplest 
polarization experiments holds in virtue of the invariance of the interaction 


+ Note added in proof. In later equations H has been taken zero since the term proportional 
to H is linearly related to the terms D, E, J. This has been pointed out to me by Drs. Ashkin 


and Wolfenstein. 
PROC. PHYS. SOC. LXV, Sac 12 
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potential under time-reversal. For collisions of identical particles, ce” must be 
symmetrical for the interchange (1, 2), which requires at once that B=C, D=E, 
F=G in (5), whence V =0 in (8) and (9). 

The complete elucidation of the scattering operator (8) requires the study 
of more complicated collisions, for instance of polarized nucleons with polarized 
target, somewhat beyond available techniques at present. Bhatia (1950) has, 
for some simple cases, studied theoretically the change in polarization of a fast 
nucleon in a collision with unpolarized nucleons, an experiment which provides 
the additional information given by Trace (6,,0%0,,0%'). 

The only term of (8) which couples singlet and triplet states is that proportional 
to V. In the collision of unpolarized nucleons, because of this term the 
polarization produced in the recoil nucleons is not equal to the polarization 
produced in the deflected nucleons, since Trace (o,0%e/") — Trace (o,0%o%") is 
not zero but is proportional to V. However, for the nucleon potentials currently 
proposed, c% is symmetrical for the interchange (1, 2), whence the two 
polarizations are predicted to be equal. An experimental test of this consequence 
would be of considerable interest. 
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ABSTRACF. ‘The polarizability of the hydrogen molecule is calculated theoretically. 
The energy of the molecule, under an external electrostatic field, is calculated by the 
variation method, using 11+10 terms and 11+9 terms James—Coolidge type wave 
functions, for the polarizability parallel to and perpendicular to the molecular axis 
respectively. By expanding the energy in powers of the field strength, the polarizability 
is determined. By using a Morse function for the adiabatic potential of the molecule, the 
0-0 and 0-1 matrix elements of the polarizability between vibrational states are calculated, 
with the result : a =7:89, Yo=2°78, o,=1:39, yo=0-90 for ordinary hydrogen lay. 
and =7°75, Yoo=2°68, %=1:13, yo:=0-71 for heavy hydrogen D, (in units of 
10~*° cm*), where « is the mean polarizability and y the anisotropy, 
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§1. INTRODUCTION 

HE polarizability is one of the fundamental quantities characteristic of 
| each molecule, and its value is often required when we want to study 
various phenomena concerning the molecule. Experimentally, the 
measurement of the dielectric constant and of the refractive index gives us fairly 
accurate values of the mean polarizability, at least in the vibrational ground state, 
but the anisotropy of the polarizability tensor is obtainable only with low accuracy. 
Especially in the case of non-polar molecules, we have no experimental means 
of determining the sign of the anisotropy. A purely theoretical calculation of 
the polarizability is, therefore, very desirable. In general this is also a difficult 
task, but in the case of the hydrogen molecule we can hope to obtain reliable 
results, since very accurate wave functions for the free hydrogen molecule, 
obtained by James and Coolidge, are fortunately available. The extension of 
their method to the case in which an applied electrostatic field is present enables 
us to find the polarizability. The results of our calculation along this line are 

reported in the present paper. 


§2. DEFINITIONS. COMPARISON OF RESULTS HITHERTO OBTAINED 

The polarizability of a molecule is a tensor, with three principal values «,,, 
%, and «,, in general. In the case of an axially symmetric molecule such as 
hydrogen, «,, is equal to «,,, if we take the axis of the molecule as z axis, so that 
we have only to know two independent quantities «,,, and «,,. We denote these 
quantities as~, anda, respectively. Itis convenient to introduce further « and y by 

Poe Lee oe yay ae er eee (1) 
x is called the mean polarizability and y is the measure of the anisotropy. 

Results of the quantum-mechanical calculation of the polarizability made 
hitherto by various authors are shown in table 1. Since all previous calculations 
were made using rather simple wave functions such as Heitler-London’s, Wang’s, 
Rosen’s, and the m.o. of Coulson, the accuracy of those results might not be 
very high. Values obtained in this paper are also shown in this table. 

Most of the authors calculated the polarizability for one or several fixed 
values of the internuclear distance R, but if we take the anharmonicity of 
vibration and the non-linearity of « in dependence on R into account, we find 
that the quantities which have direct relation to experiment are not the values of 
the polarizability or of its derivative at a fixed value of R, but the diagonal or non- 
diagonal matrix elements of the polarizability between vibrational states. «,y 
and y,, in the table are these matrix elements. It is to be noticed that there are 
fairly large differences between a, yoo and («)p, (y)z, For the sake of 
comparison we have given in table 1 also the values of a, and yo, multiplied by 
2n(vM/h)'”, which should be equal to 0«/AR, dy/OR if there were no anharmonicity 


and no non-linearity. 
§3. GENERAL THEORY 


The Hamiltonian of the hydrogen molecule in the presence of a uniform 
external electrostatic field F is given by 


1 
HoH GH, HH oR tr), eee (2) 


where H® is the Hamiltonian of the free hydrogen molecule and rj, r, are the 
position vectors of the two electrons. In this paper atomic units are used 
throughout the calculation, 
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Table 1. Theoretical Values of the Polarizability of the Hydrogen Molecule 


ea aT a DA oe ou = a, References 
0:85 0:61 0:77. —0-24 
0-52 0:77 0-60 0:25 
0:74 0°75 0:74 0:01 0:95 1:10 1-00 0:15 (3 a) 
0:67 0:71 0:69 0:04 0-80 1-08 0-89 0-28 (3 b) 
0-73 0:75 0:73 0-02 0:98 1:02 0-99 0-04 (3 c) 
0:66 0:71 0-68 0:05 0-89 1-02 0:93 0-13 (3 d) 
0-59 0-72 0-63 0-13 (4) 
0:70 0:92 (OE7/7/ 0:22 (5) 
0:66 0-68 0:67 0:02 (6 a) 
0-64 0:71 0-66 0:07 (6 b) 
0-62 0-71 0°65 0:09 (6 c) 
0-66 0-83 0:72 0:17 tesa 1-86 1:49 0°55 (7) 
0:659 0-906 0-741 0:247 T1352 1-969 1-411 0-838 (R=1-40) © 
(a ‘Loo (411)oo Ooo Yoo Cla ‘or C(o%))o1 Coty Cyor 2 
0:697 0:975 0-789 0-278 12255 2:285 1:596 1-033 for Hy = 
0-686 0-954 0:775 0:268 1-218 2-190 1-542 0:971 for De g 


C=2n(vM/h)"!*, v=frequency of vibration, M=mass of H or D. 


(1) Mrowka (1932); H—-L, perturbation. (2) Steensholt (1935); Wang, variation, one 
parameter. (3) Hirschfelder (1935); a Rosen, variation, two parameters, b Rosen, 
variation, one parameter, c Wang, variation, two parameters, d Wang, variation, one 
parameter. (4) Easthope (1936); H—L, perturbation. (5) Mulliken and Rieke (1941); 
H-L, perturbation. (6) Steensholt (1947); am.o. of Coulson, variation, two parameters, 
b m.o. of Coulson, variation, three parameters, c m.o. of Coulson, variation, five 
parameters. (7) Bell and Long (1950); m.o., variation, one parameter. In this table, the 
unit of figures in the first four columns («,, «,, « and y) is 10-*4 cm’, and that of figures 
in the last four columns (00,/0R, 00/0R, 6«/0R and dy/@R) is 10-16 cm?. 


Now the electronic ground state of the free hydrogen molecule has the 
symmetry 'X;, but when the electric field is applied the charge cloud is deformed, 
and the electronic eigen-function loses the character of &; more or less. If the 
field is not very strong, the wave function can be expressed as 


¥=0+9, 


where © is of symmetry 4, and ® is of symmetry %; or II, according as the 
field is applied parallel or perpendicular to the molecular axis. © is at least of 
the first order in F. Normalization of (3) gives 


[ @dr+ [ Odr=1, 


since {@@®dr=0. If we calculate the average value of the total Hamiltonian 
with (3), we have 

E=|OH°Odr+2F[OH'@dr+ [OH dr, (5) 
dropping terms which vanish on account of symmetry. We try to minimize 
(5) under the auxiliary condition (4), and express the result in powers of F as 


follows: B= Pode Pe ee eee (6) 


from which we can find an approximate value of the polarizability «, 
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For this purpose we express © and © as linear combinations of functions of 
the respective symmetry : 


2) a u6,0;, O a 7) iP , ee ER (7) 
i k 
The equations (4), (5) are written as 
SS C,CpS;5 + 3 bb, T x, —- ik. eo lekelele (8) 
ii kh 
E=X¢¢,}+2F Deb.Hji+ 25,b,He,  - 000: (9) 
45° i,k k,h 


where S,; = 6,0;dr, T;.,=J¢1¢6,d7. Differentiating (9) with respect to c, and 6,, 
with the use of the well-known method of Langrangian multipliers, we get 


X(HY— ES,)e+FIHgb,=0, sees. (10) 
j k 


DCHG a Poy Fear =) eens Meer (11) 


h 
In the limit of the vanishing field #0, (10) gives 


2(H}, — EO S,,)c =0. Oc Alp & (12) 


This is nothing but the variational treatment of the free hydrogen molecule, which 
has been solved by James and Coolidge if we identify 0, as their functions. Now, 
in order to find E to the second order in F, we have to find 4, only to the first 
order in F, so that we can put 6,=6,') F. This gives a system of simultaneous 
inhomogeneous equations for 5,1): 


(Ay, — LOT Jb, = — 2A, ara Lo) 
and the familiar eS of the perturbation anes gives 
B= EO). i? a (RLU SEAL Sa apts ale EOD oye (14) 
By comparing this with (6), we pally obtain the following formula for the 
TE a= —2EEG Hb. anes (15) 
The corresponding wave function is given by 
Erma Og PDO bey 8 pee caine (16) 


§4. METHOD OF CALCULATION AND THE NUMERICAL RESULTS 


As © we have adopted the eleven terms wave function of James and Coolidge 
(1933), that is, we have expressed © as a linear combination of @[00000], 
A[00020], [00110], A[10000], A[00001], 6[10200], A[10020], 410110], 
A[20000], 6[00111], [00002], where 


A i} - ry SS ee 
a[ M,N, J, K, P]= Tre gata Gl perc ap Sets ways nop" pespi = (E+ &2)}; 
Reais in (17) 


é, 7 being elliptic coordinates of the electrons and p the inter-electronic distance 
divided by the internuclear distance R; ¢ is fixed to be 0-75 throughout our 
calculation. The ¢, for this value of ¢ are given in the paper by James and 
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Coolidge, for R (atomic units)=1-2, 1-3, 1-4 and 1-5, and E® have also been 
obtained by them (private communication) as in table 2. 

For the expansion of ® we have used similar functions chosen as follows. 
In the case of «,, 4[0010], ¢[1010], g[0110], [1110], 4[2010], 4[0021], 
[0030], o[2 00 1], [1012], [0112] have been adopted, where each set of 
integers has the same meaning as M NJ K of 6[M NJ K P] (in all of our additional 
terms the quantity corresponding to P is taken to be zero, to avoid too involved 
calculation, so that they do not depend on p). In the case of «,, 


lm, nj, B= 5 CEsMéany na (Es — 1)(1 — 12) P* cos py 


+ &1"Eo""nq "no! {(E2? — 1)(1 — 2°) P” cos dalexp {— (Er +E2)}, +++ (18) 
where ¢, and ¢, are the azimuthal angles of the two electrons. Here we adopted 
410000], ¢[1000], d[0100], ¢[0020], d[0002], d[001 1], d[1020], g[1002], 
#12000]. 

The details of the numerical calculation are shown in Appendix 1. 

The results of the calculation are shown in table 3. In this table the results 
are also shown for those cases in which only three, five and seven additional terms 
are adopted, in order to see the convergency of our calculation. ‘These terms 
have been taken in the order written above, except in the 7 terms case for «,. 


Table 2 
R (a.v.) 1:2 1-3 1-4 1:5 
HE) (a.v.) — 2-3260 —2+3414 —2-3460 —2-3422 
Table 3(a). a (a.v.) 

R (A.v.) ite) 1:3 1-4 1°5 
‘Three terms 4-329 4-900 5-487 6-110 
Five terms 4-800 5-342 5-952 6-630 
Seven terms* 4-860 5-420 6:051 6:753 
Ten terms 4:876 5-451 6-107 6:862 


* Out of ten terms [0030], [0021] and 4[2001] are omitted. 


Table 3 (6). «, (a.U.) 
Ike (GAG) 12 13 1-4 bs 
Three terms 3:596 3°881 4-209 4-584 
Five terms 3-688 4-003 4°364 4778 
Seven terms 3°738 4:066 4-442 4-872 
Nine terms 3-750 4-067 4-443 4:874 
If we express our results in the form 
a=a+ b(R—R,)+c(R—R,)?+ d(R—R,)°, Rae a eee (19) 


we obtain a, 6, cand das in table 4. a and b are (x), and (0x/OR)p, respectively, 
and have been given in table 1, for the sake of comparison with results obtained 
by other authors. 


Table 4 
a b c a 
# 6:1144 7:0338 4-948 2°98 
ul hes g.s. - 0-90593 x 10-24 1:9695 x 10-16 2°619 x 10-8 2°98 
4 J (AU.) 444678 4-0424 2-803 —0:34 
1\ (c.g.s.) 0:65885 x 10-*4 1°1319 x 10-16 1-483 x 10-8 —0:31 
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The adiabatic potential of the hydrogen molecule can be fairly well expressed 
by a Morse function, at least in the neighbourhood of R=R,. According to 
Hirschfelder (1950), the numerical constants have been determined as follows: 


U=109-46(—2e-*+e2) kcal/mol, eae (20) 


x = 1-0298(R/ay — 1-401), ay) =0-529151 x 10-8cm. The vibrational wave functions 
%, for this potential can be obtained analytically (v is the vibrational quantum 
number). ‘They are slightly different for H,, HD and D, due to the difference 
of their reduced masses. After a short calculation, we can also show that 


Dey = { ti WCRER Mie dRY © wet tee (21) 


can be expressed as in Appendix 2. The numerical values of these integrals 
for the three types of hydrogen molecule are shown in table 5, and using these 
results and formula (19), the matrix elements of the polarizability are obtained 
as given in table 6. 


Table.5;, Dz(4.0:) 


H, HD Ds 

D®, 1 1 1 

Di}, 0-0427 87 0-0369 57 0-0300 81 
D2, 0-0301 72 0-0257 69 0-0206 94 
D3, 0-0045 43 0-0033 69 0-0022 16 
Vay 0 a) 0 

Di, 016703 015517 013991 
D3, 0-02419 0-01937 0-01418 
D3, 0-01727 0-01347 0-00956 


Table 6. Matrix Elements of the Polarizability 


H, HD D, 
(%)o0 65782 (0:9746) 6:5119 (0-9651) 6-4350 (0-9537) 
(a1 )oo 47029 (0-6968) 4:6674 (0:6917) 46257 (0-6856) 
Loo 5-3280 (0-7894) 5-2822 (0-7829) 5+2288 (0-7749) 
Yoo 1:8753 (0-2778) 1-8445 (0-2734) 1:8093 (0-2682) 
(4)o1 1-3461 (0-1994) 1:2275 (0-1819) 1-0828 (0-1605) 
(oL,)or 0-7377 (0-1093) 0-6774 (0-1004) 06024 (0-0893) 
Oe, 0-9405 (0-1393) 08608 (0-1276) 0-7625 (0:1130) 
Yon 0-6084 (0-0901) 0-5501 (0-0815) 0-4804 (00712) 
ae 00002666 00002643 0)-0002616 
p, (Ray.)  0:0162 0-0160 0-0157 
p, (Raman) 0-052 0-0512 0-0500 


In this table, the values of the polarizability in atomic units are given first, and their 
values in units of 10-24 cm® are given in brackets. 


We see in table 3 that the calculated values of the polarizability increase with 
the number of added terms in the wave function. ‘This is quite natural since we 
have made use of the variational method, and the energy is expressed as (14) 
in which E® is constant throughout our calculation. If the James—Coolidge 
value of E is close enough to the actual value, which seems to be true, our 
theoretical values of the polarizability are a little lower than the true ones. ‘There 
may be some additional errors in calculating the matrix elements of the polariz- 
ability, since we have approximated to the adiabatic potential by a Morse 
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function, and neglected terms higher than (R—R,)* in the expansion (19), but 
our estimation of errors from various sources shows that the final errors in («,)oo 
and (a,)o9 do not amount to more than a few per cent. The percentages of the 
possible errors in the values of y and «—%p may be somewhat higher. 


§5. COMPARISON WITH EXPERIMENT 


The dielectric constant ¢ and the refractivity n (extrapolated to zero frequency) 
of gases are given by the formula 


e—1l=n?-1=47 Na), asks ucitse (22) 


where N is the number of molecules in a unit volume. Our theoretical value of 
e or n? for the ordinary hydrogen H, at s.t.p. is 1:0002666, while the experimental 
values are 1:000273 (Tangl 1908), 1-000263 (Fritts 1924), 1-000259 (Braunmiihl 
1927), 1:000265 (Zahn 1924) for «, and 1-000273 (Van Vleck 1932) for x”. 

Bell (1942) has analysed the values of the refractive index of D, measured 
by Larsén (1936), and has obtained the result: (ag,—«p,)/%q,=13%, for which 
our theoretical value is 18-6%. 

The depolarization factor in Rayleigh scattering of natural light is expressed as 


Pn=Syqo"/(45c%997+ 77 007)- eee (23) 


Our theoretical value for this quantity, for H,, is 0-0162, while the experimental 
values are 0-022 (Cabannes and Granier 1928), 0-0274 (Rao 1927), 0-0257 
(Partharasathy 1932) and 0-009 (Volkman 1935). 

The same quantity for the vibrational Raman line can be expressed by the same 
formula if we put %, and yp, instead of % 9 and yop respectively. The theoretical 
value is 0-052, while the experimental one is 0-13 (Bhagavantam 1932). 

The ratio of the intensity of the vibrational Raman line to that of the Rayleigh 
une i vR* Ho17 + 139 42/45 

Moh fey. 


where vz and vp are the’ frequencies of Raman and Rayleigh lines respectively. 
In the case of the Stokes line due to 4046 A of Hg, the theoretical value is 0-0154, 
while the experimental one is 0-0042 (Bhagavantam 1932), which is rather difficult 
to measure accurately. 

One of the authors (Mizushima 1950) has estimated a, to be about 
(0-18 x 10-*4cm* from the intensity of pressure absorption, which agrees with 
our result. 


ACKNOWLEDGMENTS 


We wish to express our sincere thanks to the members of the Kotani Laboratory 
of Tokyo University for their discussions and financial aid, to Professor James 
for his kindness in sending his unpublished data, and to Professor Coulson for 
his kindness in communicating our paper to the Physical Society. We are also 
much indebted to Misses S. Nakadate, H. Nishikubo, M. Iwano, K. Nagai, 
A. Uno, N. Fujita and 'T. Koide for the numerical calculations. 


On the Polarizability of the Hydrogen Molecule 185 


APPENDIX! 


Most of the integrals can be expressed in terms of the James—Coolidge 
X(m,n,j,k; p), but Hy, and T,,,, can be written also with the use of another set 
of auxiliary functions. For example Hj, can be reduced as follows: 
(i) Electric field parallel to the molecular axis 
d{mnjk](z13+2)6[M NJ K P| 

= 8(R/2)"{U(m+ M+1,n+N,j+F+1,k+K; P) 
+ U(m+ M,n+N+1,j+J,kR+K+1; P) 
+ U(m+N+1,n+M,J+K+1,k+J; P) 
+ U(m+N,n+M+1,j7+K,k+J+1; P)} 

=2(R/2)'{X(m+ M+1,n+N+2,j+7+1,k+K; P) 
—A(m talon Ni jd lke +2; P) 
+X(m+ M,n+N+3,j+J,kR+K+4+1; P) 
—X(m+ M,n+N+1,7+J,k+K+3; P) 
+X(m+N+1,n+M+2,j+K+1,k+J; P) 
—X(m+N+1,n+M,j+K+1,k+J4+2; P) 
+X(m+ N,n+M+3,j)+K,k+J+1; P) 
—X(m+N,n+M+4+1,j+K,kR+J+3;P)h ae (B) 

(ii) Electric field perpendicular to the molecular axis 
d(mnj k(x, +x2.)0[/M NJ K P| 

=4(R/2)"{U(m+ M,n+N,j+J,k+K; P)+U(m+N,n+M,j+K,k+J; P)} 
+2(R/2)"{V(m+ M,n+N,j+J,k+K; P) 
+V(m+N,n+M,j+K,k+J; P)} 

=2(R/2){X(m+ M,n+N+2,j+J,k+K; P) 
—~X(m+M,n+N,j+J,k+K+2; P)+X(m+N,n+M+2,j+K,k+J; P) 
—~X(m+N,n+M,j+K,k+J+2; P)+X(m+ M,n+N+2,j+J,k+K; P) 
—X(m+ M,n+N,j+J,k+K+2; P)+X(m+N,n+M+2,j+K,k+J; P) 
—X(m+N,n+M+2,j+K,k+J+2;P)},  } }¢£ saceee (B) 


where U(m,n,j,k; P) = ; oe {I exp { —2c(E, + &5) }E:E0"Ny'2" 


x(E" — nr )(Es" a 3)p” dé, dé, dy, dns dd, dh», 
U(m, n,j,k; P) =U(m+2,n,j,k; P)+ U(m,n,j+2,k; P) 
— U(m+2,n,j+2,k; P)—U(m,n,j,k; P), 
8s 1 aMay Wee 
V(mmjk: P) = 5 a [J [Jf exp {—2elér + &a)eamearn’ne® 
x (€,2 —1)(1 — 4”) cos 264(&? — 3")(Eo® — 2”)p" dE dE dy, dy dd, ddy. 
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U, Uand V can be integrated elementarily in those cases when P=0 or 2. Inthe 
case when P=1, the Neumann expansion of 1/p can be used, and, for example, 
U(m, n,j,k; 1) is expressed as follows: 


U(m,n,j, ks 1)= E D,W,(m+4,n+2)Cj)C(k) + W,(m+2, n+4)C (CR) 
—W,(m+ 4, n)C(j)C(k+2)—2W, (m+ 2,n+2)C,(j)C(R) 
—W,(m,n+4)C(j+2)C(k) + W(m+2,n)C(7j){2C(R+2)—ClR+4)$ 
+ W,(m, n+ 2){2C,(j +2) — C7 + 4)3C(R) 

— W{m, n){2C,(j+2)C(k+2)— C+ 4)C(k+2)—CG+2)C,(R + 4)} 
—2W,(m+3,n+3)C(j+1)C(k+1)+2W(m+3,n+1)C(7+1)C(k+3) 
ra es n+3)C(j+3)C(k+1)-2W(m+1,n+1)C(j+3)C(k+3)] 
- 2 DUW2(m+2,n+2)C1(7)C MR) — W3(m + 2, n)CA(J)C (Rk +2) 
_ v3 Vm, n+2)C3(7+2)C3(k) + W(m, n)C3(j + 1)C A(R +2), 

in which D,°=27+1, D3 = —2(27 + 1)/7?(7 +. 1)?, and 


C20) =4 [ PGa)(1—n) ado, 


W,*(m,n) =| | exp {—2clEs + &5)}0,"Ex)P,(E ME? — (Ea! — I EMEa" dE, dl, 


€,(€_) being the larger (smaller) one of €, and €,. Superscript v is omitted if it 
is zero. 
We carried out the calculation of Hj, and T,, by both methods mentioned 
above, and, by comparing the two results, we could avoid numerical errors. 
The numerical values of the integrals required in our calculation have been 
partly published by Kotani, Amemiya and Simose (1938) and Kotani and 
Amemiya (1940), and the remaining part may be published elsewhere. 


APPENDIX 2 
The vibrational wave function corresponding to the Morse function 
U=D(—2e-*+e%*), x =a(R—R,) 
is given by 


a v ac, OB ae 
— J gi e 2 g(k-)2 y Sot AES a Deg 
by v s=oS8!(v—s)! I(k—v-—s)’ 


where z=kexp{—a(R—R,)}, k=4n(2uD)1"/ah. With this wave function we 
can show that 


a+ 0 
ie | thy(R — R,)”by aR. 
This can be reduced to the following formulae: 
Dy =1, Dy = (vAv') 
1 1 (A= 3)e2 5 
LV ae a ne ice oe 
Doo i {logs WO) + gq . Do = a | a Ys 


Dio Za (louk HA) + p37) + gap tH 


——— 
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(k—3)2 2 7 1 1 
imetepee =a logk YA) + pz + gap. 


D3, = -a{(logk—w(e) 4 mi) ay, (Iogk— Hh) + ri} 


2 
+ ap +3 (lok H+ 5) HH) yay}, 


Sea mats Ree Ise Lo? 
$s Saar aera logk—Ylh) + gy + Fy) 


where logk-W(8)= 55 +55 — gn toe 
Pa Ue tee 
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ABSTRACT. It is shown that on the basis of the Bloch—Peierls—Bardeen theory 
of conductivity, the use of 0g—@p in the Bloch—Griineisen formula is a reasonable 
approximation in the temperature range TZ Op. At very low temperatures T < Op, 
however, 6g=6;,. This, it is shown, disagrees with experimental data much more than 
the ad hoc use of 0g =O@p in the Bloch—Griineisen formula. 


§1. INTRODUCTION 


s far as we know, in all theoretical and experimental considerations on the 
electrical conductivities of metals, it is tacitly assumed that @ (in future 
denoted by @,) is equal to 6p, the Debye characteristic temperature. 

Although this assumption has received remarkable experimental support—at 
least in the temperature range 7 ~0,—mainly from the work of Griineisen 
(1933), so far its justification on theoretical grounds has not been attempted. 
Indeed, according to Bloch’s (1928) treatment, where only the longitudinal 
elastic waves contribute to electrical resistance, one would expect 0g=0 
corresponding to longitudinal elastic waves (which we shall denote by 6,).T 
However, Blackman (1951) has shown that values for 6; for most metals are 
appreciably larger than those for 0, or 0». 

It is the purpose of this paper to show that the conclusion 6, = 0;, is considerably 
modified when one also takes into account the so-called Umklapprozesse, which 
Bloch’s treatment neglects. In fact, we shall see that the use of 6,6) in the 
Bloch-Griineisen formula is a fair enough approximation for T2@p. For low 
temperatures (7’<@,), however, one should put 0, =0;. As we shall see in §3, 
the use of 0,=06; in the Bloch—Griineisen formula gives much more violent 
disagreement with the experimental data than the ad hoc use of 0,0). It is, 
therefore, felt that the assumptions of the theory for low temperatures need 
re-examination. 

It may be mentioned here that, since any exact treatment of these topics is 
forbiddingly difficult—mainly because of the intractability of Bloch’s integral 
equation under general conditions—our considerations are essentially of a 


semi-quantitative nature. Moreover, this discussion will be confined to 
monovalent metals only. 


* Now at the Department of Mathematical Physics, University of Edinburgh. 


+ We wish to mention here that Dr. Blackman’s (1951) emphasis on this point is, to a large 
extent, responsible for the present paper. 
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§2. THE CASE T2 4 

In this section we consider the case T2605. However, before proceeding 
further, we mention that the expression for the electrical resistance o is derived 
in three steps, (a) the calculation of the probability P(kk’) dS(k’) per unit time 
for an electron in state k to jump into state k’ in an area dS(k’) of the Fermi 
surface, (b) to solve Bloch’s integral equation for the mean free path J, and 
(c) to obtain o, which is simple once / is known. 

First we consider explicitly the case T>6). Then P(kk’) is given (cf. 
Bardeen 1937, Krishnan and Bhatia 1948) by 


3 2 
P(kk’) = F(Jk—k’|)S > ee = F($)S, say, 
m j=1 ai 

in which F(| k—k’ |) is a function of k (the electron wave number at the Fermi 
surface) and the angle of scattering ¢ only in the combination ksin}¢, the 
precise functional dependence of F on ksin}¢ being determined by the 
perturbing potential, q is the wave vector of the elastic wave, n,, is the unit 
displacement vector of that wave of which the frequency is v,,, and B,, is a vector 
in the reciprocal lattice, corresponding to some lattice spacing of the crystal; 
k is related to k’ by 


k—-k’=B,+q and |k’|=[kl. 


First consider B,,=(0,0,0). If we assume the velocities V, of the elastic 
waves to be independent of q, then (1) reduces to 


cos? y; 


P(kk’) = F(¢) DS WI ECD OG 5 SAV Ae yee Ae ate (3) 


in which cosy;=(n,;.4q,;/¢;). Here then only the longitudinal component of 
the elastic waves contributes to P(kk’). In an elastically anisotropic crystal 
€=C4, —Cyg —2€4440, the factor S, in (3) depends on the direction of q in a 
complicated manner (see for example Krishnan and Bhatia 1948). Now, in order 
that the step (b) above could be taken, one assumes « =0, when 


gent ae egal nt Pe Ie. (4) 


in which q,, is the maximum possible value of g. However, in all usual treatments 
SS, in (3) is replaced by 
[3(1/Vp3 + 2/7?) P° =(hgn/*8p)? 


instead of by (4); this procedure, in view of the above, is not justified. It may 
be mentioned that when «0 one might replace S; (as an approximation) by its 
average over all directions of q but not by the average of the left-hand side of (5) 
as is usually done. 

Consider next terms for which B,,40 of (1) (Umklapprozesse). ‘These terms 
contribute to P(kk’) when ¢ is greater than about 80°. Now here, even under 
the assumption of elastic isotropy, « =0, the sum S—.S, would not be a function 
of ¢ alone. Hence any exact treatment is difficult. However, a qualitative 
estimate of the contribution from these terms to (1) can be easily obtained. 

First we note that these terms for which B,, 40 depend upon the component 
of the displacement vector along q+ B,, and not along q, as is the case for B,, =0. 
Now it is easy to see that the angles between the vectors q and q+B,, cover a 
wide range, as the directions of incidence and scattering are varied. It is, 
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therefore, reasonable to conclude that the contributions to S—S, averaged 
over all directions of incidence and of scattering (4 >80°) from longitudinal 
and transverse waves are of comparable magnitude. Hence, if the Umklapprozesse 
alone were responsible for the electrical resistance of a metal, one would expect 
@, to be very close to 65. 

Next we give some estimates for the relative contributions to the sical 
resistance of sodium metal at room temperature (7 >6,, or 6) from the B,, =0 
term and the Umklapprozesse. For a given direction of incident electrons in the 
crystal, let R(l,), [R(«)=1/x], and R(/,) denote the contributions to R(/) from 
these two terms respectively.{ Then the contributions o, and oy to o are given 


by e= 6, -oyseenstans el) Re eee (6) 


in which R(D) denotes a suitable average of R(/) over all directions of incident 
electrons in the crystal. Krishnan and Bhatia (1948) estimated the value of 
R(L)/R(h) = 09/0, =, say, by suitably integrating P(kk’) for a few directions of 
incident electrons. They obtain v~3. However, they have overestimated the 
value of v since they took F(¢) to be nearly isotropic. Using the same method, 
but taking for F(¢) the function given by Bardeen (1937), one obtains v=2. 
We have also calculated v making Bardeen’s assumptions, which make P(kk’) 
a function of ¢ alone, and have obtained v=1-2. His assumptions, however, are 
of such a nature as to lead to an underestimate of v (see Bardeen 1937, p. 695). 
Although these calculations are by no means exact, they demonstrate conclusively 
that at room temperatures Umklapprozesse contribute most to the electrical 
resistance of sodium. However, since about one-third of the total resistance 
comes from the term for which B,,=0, to which only longitudinal waves 
contribute, if one were to use a single value of 6p, it should lie intermediate 
between 6, and 0, but closer to 0, than to 6;. 

For noble metals, which are face-centred cubic and for which the ratio of the 
smallest |B,,| to the diameter of the Fermi sphere (| B,,; |/2k=1-1) is smaller 
than the corresponding ratio for the alkali metals (| By, |/2k =1-14), obviously 
the relative contribution from the Umklapprozesse would be still larger than that 
for alkali metals. Hence for noble metals the effective 6, values would be much 
closer to 0p than for alkali metals. This is demonstrated in the table. 


Metal Li Na K Cu Ag Au 

ot Gh 540 282 — 540 — —- 
On 328 150 100 315 215 170 

ip keds 363 202 163 333 223 175 


* From Blackman (1951). 
t+ As determined from resistivities in the temperature range T~Op (Griineisen 1933). 


The above arguments, qualitatively and to a large extent quantitatively, still 
remain valid at temperatures 7’~4, since these temperatures are large enough 
to excite appreciably all q-quanta. ‘Thus we see that, at least for T= 0,, 
O,~@p is a fair enough approximation. One may, therefore, take empirically 
the expression for o to be given by the Bloch—Griineisen formula: 


a/o(Ty)=ATG(Og/ Ty format cp) Winn ehac eee (7) 
where 6) <0, <0,, AT) ~1, o(T)) is the specific resistance at temperature 7) > 0, 
or 0 and G(y) is the Griineisen function of the variable y. As mentioned before, 


t Although the separation of two contributions to R(l) may not be possible under general 
conditions, we assume for convenience that this separation is possible, 
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(7) is in very good agreement with experimental data. An empirical formula 
for o more appropriate than (7) would be (7'= 0) 


o/o(T>) = AT[xG(O,/T)+(1—«)G(0p/T)], we eee (8) 


where x =0-3 to 0-4 for alkali metals and less than this for noble metals. 

One should remember here, however, that the Bloch—Griineisen formula is 
obtained under very restrictive assumptions which are not justified when one 
takes Umklapprozesse into account and (7) and (8) are more or less only empirica] 
formulae. A detailed investigation regarding this point is now being carried out, 


$3. THE CASE T<ép 


In the temperature range 7~0), (since 6)<@,) the second term in (8) 
decreases less rapidly with decrease of temperature than the first term. However, 
at still lower temperatures, the contribution to c from Umklapprozesse will not 
be proportional to the Griineisen function G(6,/7T) for the following reason. 
As is obvious from (2), the smallest elastic wave number g, which can give rise 
to Umklapprozesse is given by q,=(|B|—2k), where B is the smallest of the 
vectors B,, of the reciprocal lattice. Hence, the temperature 0, below which the 
Umklapprozesse become almost completely ineffective may be taken as 0, ~(q,/qm)4p- 
(For sodium and lithium, to take two examples, @,~30° k and 70° respectively.) 
Therefore, at temperatures below @,, the first term of (8) alone will contribute 
too. In this range of temperatures, obviously, this term would give both a very 
much sharper decrease of o with decreasing T and a smaller absolute value 
of othan(7). Recently MacDonald and Mendelssohn (1950) have made a detailed 
comparison of the observed resistivities of the alkali metals at low temperatures 
with the Bloch—Griineisen formula (7). They find that, although in the 
temperature range T~6, the agreement with (7) is good when @, is chosen in 
accordance with the table, at very low temperatures (T<0,) formula (7) gives 
for all the alkali metals (except for sodium where the discrepancy is small) a much 
sharper decrease of o with T than is actually observed. For example, for K the 
value of 0, needed in (7) for agreement at these temperatures is about 70° k which 
is even smaller than 0). In view of the above discussion, this discrepancy between 
calculated and experimental results is very much accentuated when one uses 
for o, as one should at these low temperatures, the first term of (8) alone instead 
of (7). 

Thus we see that, although the Bloch—Bardeen theory is satisfactory in the 
temperature range T'= Op, it is in serious disagreement with experiments at very 
low temperatures. It may be mentioned that at these low temperatures the 
observed specific heats show a similar deviation from the Debye formula. 
Whether the discrepancy in the Bloch—Bardeen theory is due to its use of Debye 
vibrational spectrum or due to some other cause remains to be seen. 
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ABSTRACT. A calculation is made of the energy density of electronic states for the 3d 
electrons in metallic nickel, using the approximation of tight binding. The early stages of 
the calculations and preliminary results have been previously summarized by Fletcher and 
Wohlfarth in 1951. The theoretical basis of the approximation is outlined and the validity 
of the underlying assumptions examined. The required overlap integrals are evaluated 
and the secular equation for the energy as a function of wave vector is solved exactly over the 
energy range of interest for nickel. The density of states curve is computed graphically. 
Among results of physical interest obtained are : bandwidth 2:7, ev, degeneracy temperature 
1:4; x 10® deg k, low temperature electronic heat coefficient y 1-6, x 10~* cal. mol deg 
(compared with the experimental value of 1:74%10-*). The results are compared with 
those obtained by the cellular method of approximation used by Krutter in 1935, and 


Slater in 1936. 


§1. INTRODUCTION 

HE transition metals are of particular interest in the theory of metals, 
| since they display anomalous electronic properties, which are only 
explicable on the assumption that the energy density of electronic 
states N(£) for certain critical values of energy E is much larger than for other 
metals (Mott 1936a, b). Results deduced from the only previous calculation 
of N(£) for the 3d electrons in nickel (Krutter 1935, Slater 1936) do not agree 
very well with experiment and in the present work a re-calculation has been 
made using the approximation of tight binding (T.B.) (Bloch 1928 and others). 
The early stages of the calculations and preliminary results have previously 
been summarized (Fletcher and Wohlfarth 1951, to be referred to as I), and the 
purpose of the present paper is to give details of the complete work, together 
with a critical examination of the underlying assumptions and a discussion of 
the results obtained. Approximate wave functions for an electron in the metal 
were formed from the five atomic 3d functions and the resulting fifth-order 
secular equation has been solved exactly, whereas previously it has been the 
practice to make certain arbitrary approximations for the sake of simplicity 
(e.g. I, and Jones and Mott 1937). These approximations are shown to be invalid. 
In proceeding with the solution of this equation, numerical calculations have 
been made for the fi.st time of the overlap integrals occurring in the matrix 
elements of energy. A feature of the calculation is the occurrence of the correct 
energy degeneracies as predicted by group theory. Certain results of physical 

interest are discussed in §4. 


§2. THEORY OF THE APPROXIMATION 


The complete theory of the 'T’.B. approximation has been given by a number 
of authors (e.g. Mott and Jones 1936) and will not be repeated here. Bloch 
wave functions are built up from suitable atomic functions; those used here 
correspond to the five-fold degenerate 3d level of the isolated atom, i.e. 


$1 =(15/4m)Pxyf(r)/7", y= (15/4m)¥2yaf(r)/7? 
$= (15/4r)Paxaf(r)/1", b4=(15/16m)"%(a?— yf)? (1) 
bs =(5/167)19(322 — 9 )f(r)/7? 
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where f(r) is the normalized radial function of the isolated atom and the 
numerical factors are normalizing constants for the angular functions. From (1) 
the functions 


Par => exp (zal. k)¢,(r—al), n=1....5 ale ake (2) 


(where k is the wave-vector and al the position of any atom in the crystal lattice, 
summation being over all atoms) are constructed and may be shown to be of the 
correct Bloch form exp (ék.r)u,(r), where u,(r) has the periodicity of the lattice. 
For a cubic lattice a is simply the lattice constant and J,, l,, 1, may be integral 
or half-integral. Finally, the zero-order wave function for an electron in the 
crystal is assumed to be a linear combination of the functions (2) 


¥ (Kk sr)=d Fnm(K) bm x(¥)« GO.0006 (3) 
m™ . 
If H= —(h?/2m)V?+ V(r) denotes the Hamiltonian and E the energy for the 
electron in the crystal, then Y’,, must satisfy 
(ie Ey OT al Gli ae Be lis muah ear ete (4) 
If the non-orthogonality of the %,,, with respect to n is neglected, an approximation 
to be discussed later, it is easily shown that E must satisfy the secular determinantal 


con Tere (er ts ett ee ee (5) 
where H,,,,=[ftyp* Hing 7) /On%m, and a,?2=fb,,*%,,d7r. It can readily be 
shown that 

a,2=N¥ exp(—ial.k) | $,*(r—al)$,(r)dr=N, ss (6) 
where N is the total snitare of atoms in the crystal, on the assumption that the 


contributions to (6) from all other atoms (140) is negligible compared with that 
from the central atom (1=0). Finally, H.,,,, can be expressed in the form 


Fim =(Eo+C) nm + = exp (—ial -k) | $,*(r—al)(V—U)b,(r) dr, oo... (7) 


where U is the potential due to an isolated atom at l=0 and @& is the energy of 
a 3d electron in the field of such an atom. The constant C, termed the self-energy 
of the central atom, will be defined below. For convenience, energies may be 
measured above and below the constant value 6)+C. With the assumptions 
to be discussed below, the matrix elements H.,,,, are then found to have the 
following explicit values 


H,, = —4A, cos € cosy +4A, (cos 7 cos f+ cos £ cos €) 

H,, = —4A,cos7 cos € + 4A, (cos € cos € + cos € cos n) 

Hi, = —4A,cos ¢ cos € + 4A, (cos € cos 7 + cos 7 cos £) 

H,,=4A, cos € cos 7 —4A; (cosy cos + cos £ cos €) | 

H,; = — (4/3)(A,+4A;) cos € cos n + (4/3)(2A, — As5)(cos 7 cos ¢ + cos ¢ cos €) 

H,. =H», = —4AgsinCsin€; Hyg = Hy, = —4A,sin€ siny L 

Hy, =H,,= —4A,sinnsin’€; Hy,=Hy, =0 is 

Hy, = Hy. = —4Agsingn sin’; H3,=Hyz=4Agsin ¢ sin € 

As = Hs, = — (8/31) Ag sin€ sinn; Hy; = Hs, =(4/3"?)A, sin 7 sin ¢ 

Hy5 = Hs3=(4/3"?) Ag sin Csin€ | 

Hy; = Hq =(4/3"?)(A, + A;)(cos y cos €—cos ¢ cos €) j 
Mane (8) 
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where é =4ak,,; 1 =4ak,; ¢=3ak, and 
Ay=~ [$,8(e— Jay 4a, 2(V—U) bly, 2) dr | 
Ay= [bs%(sy—30s—Ja(V—Uybssy,2)dr 
As= | $.*(v—4a,y, 2— 44)(V — U) polo 9 2) dr : 
Ay= | bs*(x— $a, y—44,2)(V— Ub y,2)dr 
As=— | da*(,9— ba, 2—44)(V — U) daly, 2) dr 
Ag= | ba*(x,y —4a, 2 3a)(V — U)h4(x, ¥, 2) dr. 


The assumptions inherent in the above treatment will now be briefly 
considered. 

(i) The fundamental assumption of ‘tight binding’ is that the atoms in the 
crystal are far enough apart for interaction between them to be small. It is 
this interaction that causes the discrete energy levels of the isolated atom to 
broaden into bands of allowed energy for the crystal, and the assumption can 
only be justified on the grounds of consistency if the widths obtained for these 
bands are comparatively small. The value of 2-7, ev obtained for the bandwidth 
in the present calculation satisfies this criterion on comparison with 
bandwidths obtained for electrons farther from the atomic nucleus and therefore 
affected to a greater extent by neighbouring atoms in the crystal. 

(ii) The non-orthogonality of the y,,, with respect to m is neglected, that is, 
the approximation is made that 


EE exp fil —-1').1} [by*(r— al’) $(r— al) dr =0(1Am). (10) 


This approximation is linked up with another mentioned above, viz. that 
#,7=N and it is convenient to consider them together. As far as the derivation 
of eqns. (5) and (7) is concerned the neglect of non-orthogonality may be 
expressed in the alternative form 


[Yor tmetr < [puctbandr (nem) sss, (11) 


or¥ ¥ expfia(l—1’) .k} | $,*(r—al’)d,,(r —al) dr 
< DX exp {ia(l—1').k} i bn*(r—al’)4,(r— al) dr 
or SS exp fia(l—1').k} | $,*(r—al’)d,,(r—al) dr 
<N+ EE exp fia(l-1/). k}| hi, (P= al ein clic: ae (12) 
The integrals involved will obviously be much greater for nearest neighbours 
than for other pairs of atoms; considering nearest neighbours only for n=1, 


m=2, for example, a calculation gave for the left- and right-hand sides of (12) 
respectively the values 


—G-028N sin € sin € and N[1 +0:08 cos € cos 7 — 0-02 (cos n cos. + cos t cos €)] 
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so that both approximations appear to be reasonable. It is stressed that, at least 
in the present problem, the actual orthogonality of the y,,, is immaterial, 
provided the inequality (12) is obeyed. 

(iii) It is assumed in deriving (7) that 


[ba®(Y(V-U)b(t)dr=C8pm caves (13) 


Now V—U is an even function of x, y, 3, x—y or x+y because of the 
symmetry of the lattice and hence, from considerations of the parity of 
$,*¢m With respect to one or other of these variables, it is easily shown that 
S¢,*(r)(V — U)¢,,(r) dr =0(nAm). Also V—U must be symmetrical in x, y and z 
so that | ¢,,*(r)(V — U)¢,(r) dr will have the same value for n=1, 2 and 3. An 
approximate calculation shows that the absolute values of this integral for 
n=4 and 5 are slightly less than that for n=1, 2 and 3 as would be expected 
since ¢,, ¢, and ¢3 have antinodes in the directions of nearest neighbours while 
¢4, ?; have antinodes in the directions of next nearest neighbours. Owing to 
the assumption (see §3) that V=U for r<a/23?, however, the differences are 
very small and the assumption expressed in (13) is unlikely to introduce any 
serious error. 

(iv) The usual assumption was made that the major contribution to the H.,,, 
arises from nearest neighbouring atoms and only these were considered. Any. 
given atom in the face-centred cubic lattice has twelve nearest neighbours given 
by 1=(+4, +4, 0), (+4, 0, +4), (0, +4, +4) and six next nearest given. by 
I=; +1, 0, 0), (0; +1, 0)F (0, 0, +1). It seems probable, ‘therefore, that: the 
latter will also make a fairly large contribution, since they are only \/2 times as 
far away as the nearest neighbours. This approximation seems the most likely 
to introduce serious error but is unavoidable without seriously increasing the 
complexity of the expressions for the H,,, and hence of the calculation. 


$35 DETAILS OF COMPUTATION 
In order to calculate the numerical values of the coefficients A; some 
assumption must be made regarding the radial wave function f(r) and the 
potential Y—U. For the latter it was assumed that (i) V =U within a sphere 
about the central atom of radius A equal to half the distance between nearest 
neighbours (A =a/2**) (Mott and Jones 1936, p. 67), and (ii) U=0 elsewhere, 
For V outside this sphere an atomic potential centred on the nearest neighbour 
was used, which should be a good approximation in the regions where the 
product of the two wave functions is large. Since the atomic wave functions 
and potential have not been calculated for Ni, Hartree and Hartree’s figures for 
Cu*+ with exchange (1936) were used, as this ion contains the same number of 
electrons as the Ni atom. For ease of computation as close a fit as possible was 
made to these values by explicit algebraic expressions. With primed symbols 
referring to the nearest neighbour, the range r’<2A/5 was found to be relatively 
unimportant and the expressions 
We) wee Sone exp e- or Je Lover 2expiGa 27) 0 aces (14) 
Vir y=—{1+28 exp(—3r) Hr esa (15) 


were obtained (V is in atomic units). Integration was carried out over all 
space T less the atomic sphere 7 of the central atom. Within this range 


gy eo Oamexp: (ce Cr) ee eo. ee (16) 
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is a sufficiently good approximation, while for integration over 7 it was also 
assumed that f(r’)=1:9797’2 exp(—2r’) and V(r')=—1/r’. Although the 
r? exp (—5r) term would not be small in integrating over 7, it must be omitted 
as it has been neglected in integrating over J. Thus for instance 


Ayo i [1-979 exp (—2r) {1 +28 exp (—37’)} 
x {85-88 exp (—5r’) + 1-979 exp (— 27’) }xyx'y’/r’] dT 
— [ [1-979 exp (—2r) x 1:979 exp (—2r')xyx'y'/r']dr.  ...... (17) 


In carrying out the integrations, prolate spheroidal coordinates with the 
central atom and nearest neighbour as foci were used. The values obtained 
for the A; may then be put in the convenient form 


A, =0-1928E,, A, =0:0572E,, A, =0-0776Ey, Ay =0°1348E,, 
Az =0:0247 Eo Ag = O-U862E ar +) Sed). 1d ace Pee (18) 
where Ey =4(4, + A.) = TS49 a ee eee (19) 


A qualitative consideration of the various approximations indicated above 
shows that their effects on the A, are all likely to be small, some being additive 
and some subtractive. Their total effect appears negligible compared with 
that arising from neglect of next-nearest neighbours. 

Owing to the difficulty in solving a fifth order equation with coefficients as 
complex as (5), it was decided as a first approximation to neglect arbitrarily 
jim for n=1, 2, 3; m=4, 5, as has been done by previous workers (e.g. Jones 
and Mott 1937), and solve the resulting cubic and quadratic equations. Details 
of this work have already been given in I. Apart from giving correctly the total 
width of the band (2-7,ev) which is unaffected by the neglected terms, this 
preliminary calculation served to reduce greatly the work involved in solving 
the complete quintic equation and also demonstrated the effect of neglecting the 
above terms (see § 4). 

From a study of the N(£) curve obtained in I it was decided to calculate the 
more exact curve over the range of energy values 0-5<H/E,<0-771 
(2:35 <E<2-70ev for diagram in I), a complete calculation being excessively 
laborious and unnecessary for most purposes. In effect only those solutions 
of (5) were found, which fell within this range; further reference to the previous 
calculation indicated in what regions of the Brillouin zone such solutions were 
likely to be found. 

Along certain directions in k-space it is found that the secular eqn. (5) can 
be solved exactly as follows. Solutions marked * are identical with those 
obtained in I. 


(100) direction. E,,=4A,+4(A,— Aj) cos € * 
E, =-—4A,+8A, cos é % 
E, =4A,-—84A, cos & * 
E, =(4/3)[—A,—4A;+(44,—2A;) cos €] “ 
(110) direction. E, =4A, cos? €—8A, cos € * 


E,,3=4Ay, cos? + 4(A,—A,) cos + 4A, sin? € * 
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FE, and E; given by 
E? + [(44, + 16A;/3 +4A,/3)q? + (84;/3 — 16.4,/3 —8.A4,)g]E 
+ (169°/3)[(4.4,45 + A, Ay)? — (84,4; +2A,4,+4A,A, —24,A5)9 
+ (84,A,—4A,A;)] —64A,2p4/3 =0 
where p=sin €; g=cos €. 
(111) direction. E,=4[(2A,+2A,—A,) cos? €—2A,]. * 
The other four roots occur in two doubly-degenerate pairs given by 
EF? — 4[(2A, —A,+A,—2A;—A;) cos? €+ As]E 
+.16{[(A,—2A;)(2A, —A,—A,)—2A,?] cos* é 
+[A;(A,—2A,) +4A,"] cos? €-—2A4,7}=0. 
€=7;C=0 direction. E,.=+4(A,+-A,) cos n—4A, * 
E, =4A, * 
Ey 5=(4/3)[2A, — A, + 2(A,+ A;)(1 +3 cos? n)¥?]. * 

The general effect of solving the quintic equation exactly, that is the effect 
of including the elements H,,,,, neglected in I, is to separate the solutions from 
each other and to remove some of the points of degeneracy which occurred in I 
and which were accidents of that approximation. 

Returning to a consideration of eqn. (5) for a general vector k, various 
methods have been devised for solving such an equation, either in its determinantal 
or expanded form. Investigation of these methods showed, however, that 
none was suitable in the present case owing to the complexity of the H,,,, and the 
fact that the roots of (5) are not, in general, well separated. An estimate of the 
required solution for a given vector k was therefore made from the approximate 
solutions already obtained in I. This was substituted in the determinant, 
which was evaluated, a correction was applied to the solution and the process. 
repeated until a solution correct to three decimal places was obtained. The 
symmetry of eqn. (5) reduced considerably the region of the Brillouin zone to 
be investigated, and in all some 270 roots were obtained by the above method 
and checked at symmetrically equivalent points. 

For each wave vector k, eqn. (5) has five roots, and a major problem is that 
of collecting these roots into five sets, each set belonging to one Brillouin zone 
and having an associated energy band. As a first step the energy degeneracies. 
arising from symmetry considerations, reality of the Hamiltonian and accidental 
degeneracy may be investigated as indicated by Bouckaert, Smoluchowski and 
Wigner (1936) and by Herring (1937). 

Following the arguments of the above papers the following deductions 
regarding degeneracy can be made for the present case. 

(a) Due to symmetry there should be one triply degenerate level and one 
doubly degenerate level atk=0. A double degeneracy, arising from the former, 
should occur along the (100) direction and two along the (111) direction. ‘There 
is no other symmetry degeneracy apart from a double degeneracy at the 
point k=(27/a, z/a, 0). 

(6) If the Hamiltonian of a problem is real, it is possible for degeneracy, 
additional to that due to symmetry, to occur between a wave function and its 
complex conjugate. For a face-centred cubic crystal it is easily shown that no 
degeneracy of this type can occur. 
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(c) In addition to the above, various points of accidental degeneracy must 
occur in the present case. Their prediction, however, depends on a knowledge 
of the relative energies of the five wave functions ’,, (which will, in general, 
have different symmetry properties) for certain special wave vectors, e.g. 
k=(0, 0, 0), (2z/a, 0, 0), etc., so that most of these degeneracies are probably a 
result of the approximations of the calculation. They are indicated in fig. 1. 
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Fig- 1. Variation of energy with wave vector 
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(d) The normal derivative of E vanishes over the square faces of the 
Brillouin zone and along every line joining two opposite corners of a hexagonal 
face. E must have six-fold rotational Symmetry over each hexagonal face. 
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From the solutions of eqn. (5), curves showing the variation of E with k 
were drawn for each value of k,, and each value of k, from 0 to 27/a at intervals 
of z/6a. Only one other point of accidental degeneracy (see fig. 1(b)) was 
discovered, apart from those predicted. Attempts to collect the energy values 
into sets, as mentioned above, indicated quite conclusively that there is only 
one possible way of doing this. Solutions must be interchanged at all points 


Fig. 2. Contours of constant energy for certain planes of the Brillouin zone for the high energy 
set of solutions of (5) drawn at energy intervals of 0-05 Ey. (a) €=0 plane. (b) Square face 
(f£=7n). (c) Hexagonal face. (d) ¢=n/2 plane. 

of accidental degeneracy in such a way that for every value of k the highest 

energy value belongs to Brillouin zone I, the next highest to Brillouin zone II, 

and so-on. This is a perfectly valid procedure (Bouckaert, Smoluchowski and 

Wigner 1936), and the resulting discontinuities in grad E do not affect M(E) 

since on any given surface of constant energy of k-space such discontinuities 

occur only at isolated points and not over an area. 
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As in the previous calculation in I constant energy contours were drawn in 
the planes k, =v7/6a (v=0.... 12) at intervals of 0-025 E, for each of the three 
sets of energy values involved in the energy range considered. Examples of the 
(E, k) curves and the constant energy contours are shown in figs. 1 and 2. 

It follows directly from the definition of N(£) that 


f dS 
N(B)= || rca E|’ 
integrating over the surface in k-space for which the energy has the constant 
value E. To avoid numerical integration in three variables the approximation 
N(E)x<AK/AExAK was made, where AK is the volume of k-space enclosed 
between the surfaces of energies EF and E+AE, the value so obtained for N(E) 
being assumed to be that for E+4AE. This indicates the method used for 
calculating an N(E) curve for each set of energy values, the volumes AK being 
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Fig. 3. Density of states curve for nickel. EF, electronic energy in ev, measured from zero at the 
bottom of the band. N/(E), density of states, such that the total area under the (complete) 
curve corresponds to 5 states of either spin. Fermi limit for paramagnetic, ferromagnetic 
nickel as shown. --- - - Constituent N(E) curves. 


obtained by graphical integration. It was not possible to check this stage of 
the calculation, as was done in I, by measuring the areas under the N(E) curves, 
since they are incomplete, but any error introduced is unlikely to be greater 
than in the previous calculation, viz. of the order of 1%. Finally, the total N(Z) 
curve over the relevant energy range was obtained by adding the three 
constituent curves and is shown in fig. 3. 


§4. RESULTS AND DISCUSSION 

(a) The total width of the 3d band is 2-70 ev, as obtained already in I, since 
the highest and lowest energy values occur when k=(2z/a, 0, 0) for which 
Him =9 (n=1, 2, 3; m=4, 5). On comparison of fig. 3 with the curve obtained 
in I it is seen that the right-hand peak has moved to a higher energy value nearly 
at the top of the band. ‘This is readily understood on reference to fig. 1. 
Consider the (111) direction in particular, since here the effect of the previously 
neglected elements //,,,, is greatest. Fig. 1(c) then shows clearly that the result 
of including these elements is to force the higher solutions of eqn. (5) to still 
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higher energy values and the lower ones to lower values, i.e. to increase the 
number of states near the limits of the band and decrease the number near the 
middle. This suggests also, therefore, that the left-hand peak of the N(E) 
curve would occur nearer the bottom of the band than in I. The M(E) curve 
would then resemble that given by Slater (1936), but there are important 
differences between the two curves. The width of Slater’s curve is twice that 
obtained here and the height of his peak only about half that of fig. 3. The 
effect of these differences is illustrated below. 

(5) From fig. 3 the energy difference «, between the top of the band and the 
right-hand vertical line, which corresponds to 0-3 vacant states per atom, is 
0-125ev. Then by definition the degeneracy temperature for nickel is 
T =(e/ke) x 108 x 0-125 =1-4, x.108 °K. 

This compares with 1600°K and 2400°k deduced by Wohlfarth (1949) 
from the magnetic and thermal properties respectively, and about 2000°K 
deduced as above from Slater’s curve. Although the present value does not 
compare so favourably with Wohlfarth’s figures as does Slater’s value or that 
deduced in I, this is to be expected, since Wohlfarth’s calculations were based 
on the assumption of a parabolic band form, to which Slater’s N(£) curve and 
that in ] approximate far more closely than does the present curve. 

(c) The electronic property, which provides the simplest check on an N(£) 
curve as calculated here, is the electronic contribution to the specific heat. From 
Stoner (1936), the linear, low temperature, electronic heat coefficient y is given 
by y =2-835 x 10-*RN(E)), where R is the universal gas constant and N(E)p 
is the density of states for one direction of spin at the Fermi limit. For a 
paramagnetic, one can assume zero magnetization at absolute zero so that the 
lowest states are completely filled by pairs of electrons of opposite spins. In the 
case of ferromagnetic nickel and cobalt, however, Wohlfarth (1949) has shown 
that there must be almost complete magnetization at very low temperatures. 
Thus, for nickel, 5 states per atom of one spin (+) are fully occupied and only 
4-4 states per atom of the other spin (—), instead of 4-7 states of each spin. The 
only electrons which can take part in the specific heat process are those of (—) spin 
near the top of this unfilled band. The area under the curve of fig. 3 corresponds 
to 5 electronic states per atom so that the relevant value of N(£) is given by the 
left-hand vertical line, the area to the right of which corresponds to 0-6 state per 
atom, i.e. this is the Fermi limit in the unfilled (—) spin band. The resulting 
value obtained for y is 1-6, x 10-% cal mol"! deg * as compared with an observed 
value of 1:74 x 10-8, and one of 1-1 x 10-* deduced from Slater’s N(E) curve. 
A suitable value for comparison purposes would be somewhat lower than 
1-74 x 10-*, since this figure includes contributions from the s electrons and 
from d-—s transitions (Wohlfarth 1948, 1949), possibly 1-5 x10-%. The very 
close agreement obtained here, however, must be regarded as somewhat 
fortuitous, owing to the approximations made and the fact that the Fermi limit 
occurs at a steeply sloping part of the N(£) curve so that any error in its position 
would introduce a comparatively large error in y. 

It may also be noted, by comparison with I, that consideration of the 
previously neglected elements H,,,, has brought the theoretical value of y into 
closer agreement with experiment ; these elements have been neglected in previous 
calculations of a similar nature to the present one (Jones and Mott 1937). A more 
interesting point, however, is that the ‘ferromagnetic Fermi limit’ now falls to 
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the left of the band peak. If this does in fact occur, then a number of theoretical 
calculations of other electronic properties, based on a parabolic variation of N(£) 
with E at the Fermi limit, would appear to be in need of revision. 

(d) As mentioned in §3(a) a triple and a double degeneracy should occur 
between the wave functions atk=0. This is confirmed in the present calculation 
with an energy separation between the two degenerate levels of 0-88ev. In all 
applications of the cellular method, however, five-fold degeneracy is obtained 
at the centre of the Brillouin zone, a result due to the method of approximation. 

(e) Returning to a general consideration of the N(£) curve, Wohlfarth (1951) 
points out that the collective electron treatment of ferromagnetism was 
originally developed by Stoner (1938, 1939, 1948, 1951) for a ‘ parabolic’ energy 
band, for which N(E)oc(Ey—£)"?, where Ey is the upper energy limit of the 
band. Wohlfarth carries out similar calculations for a ‘rectangular’ energy 
band, for which M(E) is constant, and shows that the observed values of several 
electronic properties of nickel lie between the values calculated for parabolic 
and rectangular bands. Use of a band form intermediate between these two 
might therefore be expected to yield results in good agreement with experiment. 
The present curve may therefore be of use as a basis for a new calculation of 
suitable electronic properties for nickel. 

Another problem, to which the present work is related, is that of 
ferromagnetic anisotropy. The view has been put forward that this is due to 
spin-orbit coupling and Brooks (1940) has treated the problem theoretically on 
these lines, using ‘tight-binding’ wave functions and obtaining qualitative 
agreement with experiment for iron and nickel. In his treatment, however, the 
matrix elements H,,,, (n=1, 2, 3; m=4, 5) are neglected (asin I), whereas it has 
been shown here that this approximation is not a valid one. It is possible, 
therefore, that the results obtained by Brooks would be considerably modified 
if these elements were taken into consideration. 
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ABSTRACT. By comparison with a calibrated neutron source, the total rate of emission 
of primary neutrons in the Graphite Low Energy Experimental Pile has been found. 
The rate of emission of spontaneous fission neutrons from natural uranium is computed 
to be 59:5+3-3 neutrons/g/hr of uranium. 


$1. INTRODUCTION 


N the experiment described here, a calibrated Ra—Be neutron source was used 
] to determine the total rate of emission of primary neutrons in the Graphite 

Low Energy Experimental Pile at Harwell (the G.L.E.E.P.). From a know- 
ledge of the total amount of uranium in the G.L.E.E.P. the rate of emission of 
spontaneous fission neutrons from natural uranium was computed. 

The G.L.E.E.P. was made sub-critical by about 0-1°% (reproduction constant 
0-999) by inserting the control rods into the pile beyond the balance position. 
The neutron counting rate recorded by a boron trifluoride proportional counter, 
which was located near the boundary of the pile, was observed as a function of 
the position of a calibrated Ra—Be neutron source along the pile axis. 

It may be shown, see Appendix, that if C, is the counting rate recorded by 
the counter with the source in the pile and C) is the counting rate recorded with 
the source out of the pile 

Cy _ 1+(5/2)(ps/P) 
Cy) T+ape (p&p ” 


where S is the neutron strength of the Ra—Be source, » is the total rate of emission 
of primary neutrons by the pile, p is the probability that a fast neutron will escape 
resonance capture in the pile, « is a constant proportional to the thermal neutron 
absorption cross section of the neutron source, p is the average neutron density 
in the pile before the source is introduced, and p, is the neutron density at the 
source position before the source is introduced. 

If an absorber which has the same effective thermal neutron absorption 
cross section as the Ra—Be source be introduced, the observed counting rate 
will be C,, where C 1 

2 


Cy 1+ape/pSp 
Thus from (1) and (2) 
Cy re C2 =~ S' ps 
CS 
The experiment consisted in observing C, and C, for various values of p,/p. 
(C,—C,)/C, was then plotted as a function of p,/p, and the slope of the line 
obtained gave S/X. Since S was known absolutely, the absolute value of & 
could be obtained. 
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It should be noted that in the derivation of the foregoing formulae (see 
Appendix) it is assumed that the neutron density distribution is the same with 
the Ra—Be source present as with it absent. This assumption is true provided 
that the pile is close to criticality. 

On the other hand, the closer the pile is brought to criticality, the larger 
becomes the value of «p,2/pXp in the denominator of eqn. (1). This means that, 
with the pile very close to criticality, the effect of absorption of neutrons by the 
source tends to become of more importance than the emission of neutrons by it. 

A compromise between these two conflicting requirements had to be made. 
An analysis more detailed than that given in the Appendix was made, and showed 
that, provided the reproduction constant of the pile were closer to unity than 
0-1%, the neutron density distribution in the pile would be the same as for a 
critical pile, both with the source present and with it absent. 

The pile was therefore made sub-critical by 0-1% for the experiment. The 
experimental results obtained showed that, with this setting, the effect of absorp- 
tion by the source was small compared with the effect of emission by the source. 
That is, the variations in C, due to changing the absorber position were very much 
less than the variations in C, due to changing the source position. 


§2. EXPERIMENTAL PROCEDURE-AND RESULTS 

The boron trifluoride counter used was 3in. in diameter and 12in. long, 
and contained BF, gas at a pressure of 6 cm Hg and argon at a pressure of 6cm Hg. 
It was located near the boundary of the pile core and reflector in order to obtain 
a suitable counting rate. The electronic equipment used with the counter was 
of a standard pattern and consisted of a linear amplifier, a discriminator and two 
scalers in cascade. With the discriminator setting used, no counts were obtained 
when the high tension supply was disconnected from the counter. There was a 
small dead time (16sec) in the electronic equipment, and all counting results 
were corrected for this. 

The Ra—Be source used was the one described by Littler (1951); it had a 
neutron strength of 9-7 x 10®(+4$%)n/sec. The effective thermal neutron 
absorption cross section of this source was also known from the source calibration 
experiment (this was obtained by running the G.L.E.E.P. at high power and then 
modulating the power by oscillating the source through the pile). 

The absorber used was a brass cylinder. This had been shown by modulation 
experiments to have the same effective thermal neutron absorption cross section 
as the Ra—Be source. 

The G.L.E.E.P. was shut down for several hours before the experiment was 
started. This ensured that there would be no contribution to the primary 
sources in the pile from delayed neutron emitters, as a result of previous high 
power running. 

The neutron density distribution in the G.L.E.E.P. had been previously 
measured by the use of indium foils (Colmer and Littler 1951). From these 
measurements the values of p,/p for each source position used were computed. 
The source was placed at standard positions along the axis of the pile, these 
positions being at intervals of two pile pitches (one pitch is 7}in.). Source 
positions near the core-reflector boundary were not used, since the theory of the 
method assumes a bare pile, and the results would be erroneous if the source 
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were nearer to the boundary than two migration lengths (~50cm). After 
changing the source position (or the position of the brass absorber) no counting 
was done for fifteen minutes. This allowed the neutron density to reach a new 
steady state for each measurement. 

The results obtained are shown in the figure. The statistical accuracy of the 
counting rate ratios is within the order of 1%, and the accuracy of the values of 
p,/p is within the order of 3%. 

The slope of the straight line in the figure was determined by the method 
of least squares and was found to be 16:9. Since S=9-7x10%n/sec, we 
may deduce the value of &, the total primary neutron strength of the pile: 
xX =5-73 x 105 n/sec(+54%). 
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The G.L.E.E.P. became divergent in August 1947, and it then contained 
12 tons of uranium metal and 21 tons of uranium dioxide (Colmer and Littler 
1951). Since that time it has had more uranium added to it, to give it an excess 
reproduction constant of 0-2°%, and it now contains 3-462 x 10*kg of uranium 
element. 

Thus, assuming that all the primary neutrons in the pile are produced by 
the spontaneous fission of uranium, the spontaneous fission neutron strength is 


16:5 +0-9n/kg/sec, or 59-5 + 3-3 n/g/hr. 


§3. DISCUSSION OF RESULTS 


The figure quoted in the previous section for the number of spontaneous 
fission neutrons emitted per hour per gramme of uranium will be in error if there 
are other sources of primary neutrons in the G.L.E.E.P. ‘Two other possible 
sources of primary neutrons are those produced by cosmic rays and those 
produced by («,n) reactions. 

The rate of production of neutrons in the atmosphere by cosmic rays is too 
small to have any effect on the figure quoted in the last section. (Apart from the 
low rate of production of neutrons, only a small fraction of those produced can 
enter the pile because it is surrounded by a thick concrete shield.) An order 
of magnitude calculation on the production of neutrons from uranium by the 
capture of negative .-mesons has been made. It was assumed that the stopping 
power of uranium for y-mesons is the same as that of lead, and for the latter 
the value quoted by Koenig (1946) was used, that at sea level the flux of 
o-mesons is one per cm” per minute (Rossi 1948), and that two neutrons are 
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produced per pz-meson stopped. It was found that the neutron production from 
uranium by the capture of j.-mesons was much less than 0-1% of the production 
from spontaneous fission. 

It has been assumed that there is only a very small production of neutrons 
by the stopping of u-mesons in graphite. This assumption is based on the 
experiments of Conversi, Pancini and Piccioni (1947), who showed that nearly 
all u-mesons stopped in graphite undergo decay and do not interact with the 
carbon nucleus. 

The rate of production of neutrons by («,n) reactions in the uranium metal 
in the G.L.E.E.P. should be very small as the metal is of high purity. In the 
uranium dioxide core of the pile there may be a contribution to the primary rate 
by («, n) reactions with oxygen. Roberts (1947) quotes the neutron yield obtained 
by bombarding oxygen with polonium alpha-particles. An order of magnitude — 
calculation has been made based on these results. Allowance was made for 
the fact that we have a uranium dioxide target and not a pure oxygen target, but 
no allowance was made for the fact that the alpha-particles from uranium are of 
lower energy than those from polonium. It was found that the figure quoted 
for the number of spontaneous fission neutrons produced per hour per gramme 
of uranium may be too high by alittle under 2%. This figure represents an upper 
limit to the error, since the yield of neutrons by the («,n) reaction may change 
rapidly with alpha-particle energy, and the uranium alpha-particles are about 
1 mev less energetic than those from polonium. 

An early measurement of the number of spontaneous fission neutrons emitted 
by uranium was made by Rotblat (1941). He obtained a value of 53 + 10 n/g/hr, 
which is in agreement with the value quoted in this paper. Other measurements 
of the number of spontaneous fission neutrons produced per hour per gramme 
of uranium have been made by Scharff-Goldhaber and Klaiber (1946) and Pose 
(1943). Scharff-Goldhaber and Klaiber quote a partial decay constant for 
uranium of 7 x 10-**sec"t on the assumption that only one neutron is emitted 
per spontaneous fission. ‘This figure leads to a value of 64n/g/hr emitted from 
uranium. Pose quotes a partial half-life for the decay of uranium by spontaneous 
fission, from which a value of 65 n/g/hr emitted from uranium may be deduced. 
These authors do not quote the accuracy of their determinations, but there must 
be some uncertainty 1n Pose’s value, since he assumes a figure for the calibration 
of a neutron source. 

Whitehouse and Galbraith (1950) quote a value for the spontaneous fission 
rate in natural uranium of 24-2+41-5fissions/g/hr. If we combine their value 
with the value quoted in this paper for the rate of production of neutrons 
from spontaneous fission we have, for the number of neutrons per fission, 
v =59-5/24-2.=2-5 + 0:2 n/fission. 

It is of interest to note that this value for the number of neutrons per spon- 
taneous fission in natural uranium is the same as the value for the number of 


neutrons per thermal neutron induced fission in natural uranium, namely 
2:5 + 0-1 n/fission (Atomics 1951). 
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THEORY OF THE EXPERIMENT USING A ONE-GROUP MODEL 


The theory of the experiment given here is very similar to the three-group 
diffusion theory given by Littler (1951), and some of the results of this latter 
analysis are assumed in the present treatment. 

For a steady state in the pile we have the following diffusion equation: 


DWV*p —p/t +ps = —p/t, Jeonsealle\t) 


a 


where p, a function of position in the pile, is the density of thermal neutrons, 
k is the reproduction constant for the pile, p is the probability that a fast neutron 
will escape resonance capture in 7°8U, s is the rate of production of primary fast 
neutrons per unit volume, and D and 7 are the diffusion coefficient and mean 
lifetime of the neutrons. 

We assume that there is a small cavity in the pile, and that initially it is empty. 
Equation (A1) then describes the state of affairs everywhere in the pile except 
in the cavity itself. 

We now suppose that a neutron source is introduced into the cavity, and that 
the pile can again be in a steady state with p changed to p, and k changed to ky. 
We have, therefore, 

DV" 9, =p /T PS = Ripe OF ae (A 2) 


If we now apply Green’s theorem in the same way as in the three-group 
treatment, and assume that the pile is bare, so that p and p, will both be zero at 
the outer boundary of the pile, we have 


. jp—=/h Pt 
p(T —pS)+ps | (pr—p) V+ ——* | pp dV=0, ...... (A3) 


where p, is the value of the neutron density at the cavity, T is the rate of absorption 
of thermal neutrons by the neutron source, S is the strength of the neutron source, 
and the volume integrals extend over the whole pile. If the reproduction constant 
is not changed, then k=,, and hence, 


(EO OS\EpsVigGep)=0,  2.s.. (A4) 


where p and p, are the average neutron densities in the pile, and V’ is the volume 
of the pile. 
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If we write T=ap,, (where « is a constant proportional to the absorption 
cross section of the neutron source) and sV =, the total primary source strength 
in the pile, then we have 

Pe%Pre—PS)+PU(pi—p)=9. seen (A5) 


If we suppose that the neutron density distribution in the pile is the same when 
the neutron source is present as when it is absent, we may write 


pee 

Pe p Co 
where C, and C, are the counting rates obtained with a neutron counter.in a given 
position in the pile, before and after introducing the neutron source. 

The assumption that the neutron density distribution is not altered by inserting 

a neutron source into a sub-critical pile is only justified provided that the pile 
is very close to criticality (Wallace and LeCaine 1946). However, calculations 
of the change in the density distribution with the pile initially sub-critical by 
0:1°% showed that no serious error is incurred by making the assumption in this 
case. Furthermore, the calculations show that the density distribution with 
the pile sub-critical by 0-1°% may be taken to be the same as with the pile critical, 
without introducing any serious error. 


Using (A 5) and (A6), we have, finally, 


Cy _ 1+(S/%)(pe/p) : 7 
Ge alanis, ce i ee (A7) 


In deriving this formula we have assumed that the G.L.E.E.P. has a uniform 
core. This is not so, since in the G.L.E.E.P. there are two regions, one containing 
uranium dioxide and one containing uranium metal. The problem should 
therefore have been solved as a two-core problem. A separate calculation has 
been carried out using a one-group two-core model. It was shown that, provided 
that ps is the same in the two cores, the result is identical with (A7). However, 
s is not quite the same in the two cores since the amount of uranium element 
in a dioxide channel is not quite the same as the amount of uranium in a metal 
channel; also the values of the resonance escape probability differ by a few per 
cent in the two cores. Even in this case, however, the calculations showed that 
eqn. (A7) is correct to well within 1%. 

One other point is worthy of mention. If in eqn. (A3) we take p=p, over 
the whole pile, we have, for the change in reproduction constant on inserting 
a neutron source into the pile, 

Fao Ny Ait Dansk 5h 
Jp* dV 

This result is the one-group equivalent of that obtained by Littler (1951). 
It should be noted that it is independent of s, i.e. independent of the strength of 
primary sources in the pile. 
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ABSTRACT. An account is given of the methods used to determine the matrix elements 
of crystal field potentials with particular reference to rare earth ions. Emphasis is laid on the 
importance of Wigner coefficients in such problems and the idea of using equivalent 
angular momentum operators is developed. For convenience in applying the results 
tables of matrix elements are included. 


$1. INTRODUCTION 

N the course of some recent work on rare earth ions in crystals it has been 

necessary to evaluate the matrix elements of certain potential functions in 

representations in which the states are eigen states of the total angular 
momentum. Somewhat similar problems arise in connection with the magnetic 
properties of iron group ions in crystals, with the difference that the usual 
representations are then those in which the states are eigen states of the total 
orbital angular momentum and it is not necessary to include spin variables. 
In the case of rare earth ions the perturbations due to the crystal fields are usually 
small compared with the spin-orbit couplings, whereas in the iron group the 
inequality is reversed. This difference means that the evaluation of the relevant 
rare earth matrix elements is generally more complicated than is the case for iron 
group ions. It is the purpose of the paper to give an account of a method for 
obtaining their values in terms of one-electron radial integrals. 


§2. PERTURBATIONS AND STATES 


The potential functions which have been considered in detail are: 
V9 =X(327 —7*) 
V 2 =%(3524 — 30772" + 374) 
Vo =2(2312° — 3157224 + 1057427 — 57) 
Vf =U (x8 — 15x45" + 152x744 — y®) 
where the summations are taken over the coordinates of all electrons. It will be 
noticed that each V can be written in the form f(r) P,,"(0, ¢) and thus transforms 
under rotation according to an irreducible representation of the rotation group. 
This property is extremely useful in determining matrix elements; as a first 
application it follows that the matrix elements of V,° can be obtained from those 
of V,° purely by symmetry arguments. Alternatively, once the matrix elements 
of V,°, V,° and V,° are known, the matrix elements of all even potentials up to 
the sixth degree can be found from symmetry. This will be discussed more 
fully later in the paper. 
The states considered are formed by taking linear combinations of deter- 
minantal product states of 4f one-electron states. ‘The rare earth ions show 
PROC. PHYS. SOC.-LXV, 3—A 14 
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Russell-Saunders coupling so the product states are first combined to form states 
for which L and § are constants and these are then combined to form states for 
which J is constant. The final states can be conveniently written as |4f"; LSJJ,) 
where n denotes the number of electrons not in closed shells. The evaluation 
of the matrix elements of the potentials would be quite straightforward if the 
expressions of all states such as |4f"; LSJJ,) in terms of basic products states 
were known. This, however, is not usually the case and rather than go through 
the tedious process of obtaining them it is more profitable to determine the matrix 
elements by methods which do not require this information. 


§3. OPERATOR EQUIVALENTS 


The methods used depend largely on the result that within a manifold of 
states for which J is constant there are simple relations between the matrix 
elements of potential operators and appropriate angular momentum operators. 
For example, inside J} = constant, 


xy =oa[J,J,+I J «| 
(3524 — 30722? + 374) = B[35S,4 — 30S (J + 1)J2 +25, — 6S (J +1) + 3)7(J + 1)*], 


where « and f are numerical factors. ‘That such operator equivalents hold can 
be verified by using the Wigner coefficients to give the values of the matrix 
elements of the potential functions. The Wigner coefficients are given (van der 
Waarden 1932) by: 

neces (—1)L(j+m)!G—m)!(j'+m')!('-—m’')!(J + M)!(J — M)!}4? 
seh "G—m—v)!(g+m—A—v)!(j'+m' —v)!(7’—m'—A+v)!v!(A-vr)! 
with M=m-+m' and J=j-+)7’—A (A=0,1,2,....), using van der Waarden’s 
notation. The element <¢J,J,|V,°|J,J,) is obtained, apart from the factor « 
from c?, ,, with j =J, m=ZJ,, 7’ =4, m’ =0 and A=4, since V,,° can be regarded as 
a component of a vector in a D, space having m’=0. ‘This method can be used 
to check the equivalents once they are obtained, but it does not provide a very 
good method for determining them. ‘The Wigner coefficients will give the 
matrix elements, apart from the factor, but it is not always easy to pass from these 
to an angular momentum operator with the same elements. The following 
method is, however, available. 

Each potential transforms under rotation according to an_ irreducible 
representation of the rotation group. For an operator equivalence to hold, the 
equivalent must transform in exactly the same way, and it is a matter of determin- 
ing an angular momentum expression which has this property. There is no 
difficulty in doing this provided allowance is made for the non-commutation of 
J,,J,andJ,. For example, the two sets of functions x? — y?, 3z?—7?, xy, yz, zx 
and Lick —J 45 3S, —S(S 4 1), (J J, +S yJ q), 2S WJ, +J Jy), SS» +J ,J,) each 
transform according to D, and it is clear that they give the same transformation 
matrices. Thus x?—y? =a(J ,”—J,7), xy =4a(J J ,+J,J,) etc., within a manifold 
for which J is constant. ‘This process is readily adapted to more complicated 
expressions. Suppose an operator equivalent to «*—6x*y?+y! is required. 
This can be written as 4[(~+7y)*+(x—zy)*] and the equivalent operator is 
24l(J,, +e ,)*+(J,,—tJ,)*]. On the other hand, the operator equivalent to 
3524 — 30r*2° + 374 is not «[35J,4 — 30J,2J(J + 1) + 3J*(J + 1)?] since, for example, 
it is not correct to replace x?z* (which occurs in 773”) by J,?J,2. It should be 
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replaced by an expression consisting of all possible different combinations of 
J,,J,,J,andJ,. Some equivalences are given in table 1. 


Table 1 
Operator Equivalents inside a Manifold for which J is Constant 
X(32?—r?) =ar? [3J2—J(J +1)] 
B(3524—30r222-+ 374) = Br#[35J,4— 30) (J +1) J,2+-253,2—6) (J +1)+3J%S +1)"] 
=(23128— 31572244 105742?— 578) 
ae 231S,8—315J(J+1)J,4+ 7355.44 1053? (J +1)*J,? 
a4 Cie Rapala teat a ee tian uetieen | 


Values of «, 8 and y for Ground States of Rare Earth Ions 


Ton or B ue 
Cet*++4f1 °F, —2/35 . 2/7.45 0 
Prt4f sy, —52/11.15* —4/55.33.3 V7 1647.017-13.5.3" 
Ndt++4f* 41,,. —7 /33.33 —8.17/11.11.13.297 —17.19.5/13?.11%.33.7 
Pmt++4f4 5], 14/11.11.15 ORD Bes cations 2584 /117.13?.3.63 
Sm+++4f5 SA se 13/7.45 26/33.7.45 0 
Eut++4f8 7R, 0 0 0 
Gd+++47 8s 0 0 0 
Tbh+++4f8 ’F, —1/99 2 /11.1485 —1/13.33.2079 
Dy*++4f? "Hi 5/5 =2/9.35 —8/11.45.273 AH Nene 2aS 3277 
Hot++4£"5], —1/30.15 —1/11.2730 —5/13.33.9009 
Ert++4f11 4], 5/5 4/45.35 Piet os2 73 tah ALINE SI 7/ 
Tut++4f12 3H, 1/99 8 /3.11.1485 = 13.33.2079 
Yb+++4f13 °F), 2/63 —2/77.15 4/13.33.63 


For a single 4f electron, ~=—2/45, B=2/45.11, y=—4/9.13.33. 


§4. MULTIPLYING FACTORS 


There then remains the problem of finding the multiplying factors. This is 
done by using the fact that the potential functions do not show any dependence 
on the spin, so that similar operator equivalents hold inside manifolds in which 
L is constant. A convenient state in L, S, J,./, quantization is chosen and 
expressed in L,$, L,,.S, quantization and an equation obtained by equating the 
two expectation values of the function. Another state in L, S, L,, S, quantization 
is then chosen and expressed in one electron product states and again expectation 
values are equated. In this way sufficient relations are obtained to determine the 
factors in terms of radial integrals over 4f wave functions. ‘The technique for 
doing this is best demonstrated by means of an example. 

The ground term of Sm+++ is 4f°; °H5/ and the problem of finding the 
matrix omens of V,9= (3524 —30r*2"+ 374) inside the manifold of states 
spanning J=3 is considered. From the operator soe ase nee it follows that the 
matrix ne of V,° are equal to those of 


of 35J,4—30.$.3324257,7-6.$.3+3.(8)?.($)7]. 
The most general expression for |L=5,S=$,J=3,J,=$) in L,S,L,,S, 
§ p 
quantization 1s: 
4 d|2 ES HAL ESEF|0,B)pher= cleo! vena! (1) 
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where a? +8?+¢2+d?+e+f?=1 and all states have L=5 and S=3. Thus, 
using table 3, and noting that column F gives multiplying factors (see remark 
near end of paper) 
60a = <J= $j. §|V,°|J= $d, = 3) 

=@(L,=5,S,= —$|V,9| L,=5, S,=$) + 844, -31Vol4, — 2) +-- 

=a2{5|V9|5)+32¢4|[V.9]4)>+....+f2(0[V.°] 0), 
where, in the last line, the value of S, has been omitted since V,° does not contain 
spin variables. Inside a manifold for which L=5 

V P= p35 L,4 —30.5.6L,7 + 25L7—6.5.6+3.(5)*.(6)?] 

so that 60«=420f[6a? — 6b? —6c?—d?+4e?+ 6f?] using values from table 2 


fs =p 


appropriate to L=5. The state |L,=5,S,=$) is equivalent to {321 0-1} 
using the usual notation for determinantal product states. Hence 


++++4+ 4 r++t+ 4 
(5, 8|V.2|5,3)=25208 ={(3 21 0—117,°%3 210-1} 
= y[180 — 420 + 60 + 360 + 60] = 240y 
since inside a manifold for which 1=3 (f-electrons) 
35,1 — 30722? + 374 = y[351,4 — 30.3 .4.42+ 251? —6.3.443.(3)?.(4)7]. 
Also y.180 = (1=3, 1, =3|3524— 30772? + 374|1=3,1,=3)=27* where the wave 
function (not normalized) corresponding to |1=3,/,=3) is LE. 0) exp 374, 
B= 240 240.8 4 4 
25207 = 2520,11,180. 21.45/11 


There remains the problem of determining the coefficients a,b,c...etc. in (1). 
A convenient way of doing this is to operate on both sides of (1) with 2L.S. This 
operator has the property of commuting with J and inside a manifold in which 
L, S and J are constant it is equivalent to [J(J + 1)—S(S+1)—L(L+1)]. Inthe 
present case 2L.S= —30. ‘Then, for example, 
— 30a (L, RAS) —$|—30]{a]5, —$)+6]4, —8)+..3 

=(L,=5, S,=—§|2L.$| {a]5, —$) +5]4, 8) +...) 

= —25a+ /500. 
Enough equations can be obtained in this way to determine a, b,c... etc., the final 
values being 


P= 6 \12 a op 12 = 4 \1/2 4 1 \12 1 \1/2 1 \12 
an(Tr) > *=—(a) > om (g3) > =) eC) fag) 


These values give 
ae pee bart 6 é) eet 1 1 
60% = m145-11" 420 6. i —6. ii —6. 33 7 29 +4. 55 +6.o7 | 


- a 
Be eR 


In this way the matrix elements of V,° inside J = for Sm+++ can be obtained 
in terms of the average of r* over the radial part of the 4f wave function. A similar 
procedure is clearly available for V,° and V,°. 


Sa nits 


or 
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At this point it is convenient to make several general remarks. In the example, 

the passage from L,,.S, quantization to /,,s, quantization by 
++++ + 
| L,=5, S,=$)=3 210-1}, 
is readily made because the determinantal state is the only state for which L,=5, 
S,=$. It is a property of the ground multiplets of all ions obeying Hund’s 
Table 2 
Values of 3J,2-—J(J +1) 


Af F 
ee 44/2 43/2 45/2 £7/2 $9/2 +11/2: 413/24 15/2 
1/2 0 0) 
3/2 3 — 1 1 
52 2 — 4 — 1 5 
7/2 Se ==t5 — 3 1 7 
9/2 6 — 4 — 3 — 1 2, 6 
gi 2 1 —35 —29 —17 1 D3) 45 
13 /2 6 — 8 7 5 2 2 7. 13 
15/2 3 —21 —19 —15 — 9 — 1 9 Dal 35 
a F 
5 0 +1 +2 aia +4 +5 +6 ae i +8 
@) 0) 0 
1 1 2 1 
2 3° = "2 — 1 2 
3 3 — 4 — 3 0) 5 
~ 1 —20 —17 — 8 Wf 28 
5 3 10 9 6 — 1 6 5 
6 3 14 13 10 5 2 lal DP) 
7 1 —56 —53 — 44 —29 — 8 19 52 91 
8 3 —24 —23 —20 —15 — 8 1 12 25 40 
Table 3 
Values of 35/,4 —30J(J + 1)J 2+ 25J,2 —6J(J +1) +302(J +1)? 
Af F 
if Ose 5 12 2 972 FE 112041312 £152 
d bi 2 0 0 
3/2 0 0 0 
5/2 60 2 — 3 1 
vee wel i er ee 7 
9/2 84 18 3 —17 — 22 18 
11/2 120 28 12 18 — 33 — 27 33 
13 /2 60 108 63 —13 — 92 —132 — 77 143 
t5/2 60 189 129 23 —101 —201 —221 — 91 273 
ff F 
I= 0 1 2 3 4 5 6 7 8 
0 0 0 
1 0 0 0 
2 iy 6 — 4 1 
3 60 6 i — 7 8) 
+ 60 18 9 —11i — 21 14 
5 420 6 + —1—-—- 6 — 6 6 
6 60 84 64 11 — 54 — 96 — 66 99 
7 12 756 621 251 —249 —704 —869 —429 1001 
8 420 36 31 17 — 3 — 24 — 39 — 39 — 13 Ey? 
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rule that such a relation holds, for the rule states that the ground multipiets have 
maximum $ and maximum L. The state with maximum S, and L, is then always 
particularly simple in /,,s, quantization. In the case where no such obvious 
relation exists, such as for excited states, the passage from L,, S, to /,,s, quantiza- 
tion is more difficult and is somewhat analogous to the passage from J, to 
L,, S, already considered. It is suggested that the operator =1,.1, could be 
used for this in a process similar to that in which L.S was used. 


Table 4 


Values of 231/,6—315J(J + 1)J,4 + 735J,! + 105J(J + 1)9°J2—525)(J +1)J2 
$2947, —5I3(F +1)? + 40I2(F + 1)? -60I(J +1) 


J F 
J,= BT P=E3° “R5)P 7D CES Pete i eee 
1/2 0 0 
3/2 0 0 0 
5/2 0 0 0 0 
Fiz 1900 = 1S eas 1 
9/2. 5040 — 8 6 10 == tt 3 
1472" -7560.—. 20 4 25 112 eS 11 
13/2- 2160 —200. =25 185. 207 1 — 19, 
15/2713860.— 75. —=25 45 87 501, — 39 2 di7 65 
iPemepe 
J = 0 1 2 3 4 5 6 7 8 
0 0 0 . 
1 0 0 0 
2 0 0 0 0 
3a iSO) 296 Notes as 1 
Aes1 260, = 820 1 pee 4 
B50 40 pal 36 202-48 15 
6 7560 = 40), 690 22 43 gr naass 2 
Te 27800 2008 = 125 50; =. 497: ©, 1768 Ar=55 ve eese 4g 
8 1386012120 aae5 2 93° v1.28 65° = 78> — 160 Oe 


Having discussed the determination of matrix elements coupling states in 
the same J manifold the next step is to consider elements joining states in different 
J manifolds. It is not now possible to use operator equivalents and it is usually 
most convenient to use the Wigner coefficients directly. The variation of the 
elements with J, is readily obtained and there remains again the problem of the 
multiplying constant. 

Suppose for example that ¢J,J,|V,°lJ+2,J,) is required. Regarding V,° 
as the component of a vector in a Dg space with m’ =0 the coefficient of | J,J,) 
in V,°|J+2,J;) is required. This is equal to a constant multiplying c}_ 9 with 
j =6, m'=0, j=J +2, m=J, and A=8. The value of the constant is obtained 
by evaluating the matrix element coupling a chosen pair of states in a manner 
analogous to that already described. 

Thus the problem of reducing the matrix element of any potential function 
which transforms like the component of a vector with m’=0 in an irreducible 
representation space to a radial integration is, in principle, solved. It is often 
convenient to work with the components having m’ =0, but the reasoning is easily 
extended to potentials other than these, since once the multiplying factor has been 
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found the variation with J, is exactly that given by the corresponding Wigner 
coefficient. This dependence can be written as follows. Suppose V° and V* 
are two components of a vector in a D, space. Then 


{J,J,|V°| Vos) Here, ww with 7 =3,m' =0,j=J' ;m=J,,A=s+J =, 
and 


Ch I-+R|VE|NS,) =a, with j’=s,m'=k,j=J',m=J,A=s4+d' —J, 


where the «’s are the same and depend on Jand J’ but not onJ,. It will be noted 
that V* only has matrix elements between states for which AJ,=k. 'There is no 
difficulty in recognizing that a potential forms, say, the V" component of a vector 
in D, for, this being so, it must transform under rotation in the same way a9 
P,*(@,¢), an associated Legendre function, which is such that a rotation through 
% about Oz multiplies it by exp (zk). 

In the accompanying tables are given the operator equivalents of certain 
functions together with the corresponding matrix elements. The elements for 
a given J value have common factors which are given in the columns under F. 
The actual elements are obtained by multiplying the elements in the table: by the 
factors in the corresponding positions in column fF. In table 1 are given the 
factors for V,°, V,° and V,° appropriate to the ground states of all the rare earth 
ions. 
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The Possible Effects of «-Mesons in the Cosmic Radiation 
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ABSTRACT. The decay and nuclear interaction schemes of «-mesons are discussed. A 
tentative interpretation is made of some anomalous interactions of penetrating particles 
observed recently below ground. Further possible experiments testing the hypotheses made 
are suggested. 


§1. INTRODUCTION 


T is generally accepted that the -mesons form the penetrating component of 
| the cosmic radiation. They are thought to arise from the decay of 7-mesons, 
the latter being produced in energetic collisions of the primary radiation with 

the nuclei of nitrogen and oxygen atoms in the atmosphere. 


* On leave from Osaka City University, Japan. 


216  §. Hayakawa 


There are, however, several properties of the cosmic radiation observed below 
ground that do not fit in well with this picture. Thus, there is the question of the 
variation with depth of the intensity of the penetrating component. This has been 
discussed by Hayakawa and ‘Tomonaga (1949) on the basis of the above picture. 
These authors find it difficult to obtain a satisfactory agreement with observation 
at all depths. Also there is the existence of the pairs of penetrating particles 
discovered underground by Braddick and Hensby (1939). Observations by 
George and Trent (1949) and by Braddick, Nash and Wolfendale (1951) have 
shown that if these are produced by pu-mesons, then the cross section for their 
production is 5 x 10-?®cm? per nucleon, a figure which does not agree at all well 
with that found for the nuclear interaction of u-mesons by George and Evans (1950), 
viz. 5 x 10-8°cm? per nucleon. Finally, several cloud chamber pictures have been 
obtained underground by George and Trent showing anomalous scatter of 
single penetrating particles on passing through lead plates in the chamber. 
Assuming the particles to be 4.-mesons, the cross section for anomalous scatter is 
found to be of the order of 2 x 10-*8 cm? per nucleon. Such a large cross section is 
difficult to understand in view of the known weak interaction of .-mesons with 
nucleons, and in any case is about 100 times larger than the values found by Amaldi 
and Fidecaro (1950) at sea level. 

An investigation has been made of the extent to which the above anomalies 
might be attributed to the effects of the x-mesons recently discovered by O’Ceallaigh 
(1951). ‘The results are briefly reported below. 


§2. THE OBSERVED PROPERTIES OF k-MESONS 

According to O’Ceallaigh (1951) the following properties of «-mesons have been 
found using nuclear emulsions: (i) the mass is approximately 1 000 times that of 
the electron ; (ii) the x-meson decays at rest into a secondary charged particle of mass 
200 to 300 times that of the electron; (iii) the momenta of the secondaries are not 
unique, various values up to 250 Mev/c having been observed ; (iv) in one example, 
the secondary decay particle was brought to rest in the emulsion and decayed in 
turn into a relativistic tertiary particle; in this case the secondary particle may be 
identified as a y-meson; (v) the lifetime is not established with certainty, but 
probably lies in the range 10-!° to 10-® sec; (vi) the frequency of x-mesons observed 
to come to rest in nuclear emulsions exposed under 30cm Pb at the Jungfraujoch 
is approximately 10 °% of that of 7—y decay events. 


§3. THE INTERACTION SCHEME OF «&k-MESONS 
In the light of the above properties, the following interaction scheme of the 


k-meson is tentatively suggested. 
3.1. p-—Decay 


The y-decay of «-mesons appears to be similar to that of the B-decay of 
p-mesons : 

Keak eee | I Re eee (1) 
where both («*, «°) and (4*, 2°) behave according to pair-theory, whatever their 
values of spin. ‘The mass of the x°—particle may eventually be determined from 
the properties (1) and (iii) in §2, though at present the results are not sufficiently 
accurate. ‘The decay of the charged x«-mesons via a nucleon pair must be allowed 
as an alternative possibility, but it can be at least partially forbidden by selection 
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rules (cf. Fukuda and Hayakawa 1951). The coupling constant for the process (1) 
is of order of magnitude similar to that for the 8-decay of nucleons or y-mesons, 
a more precise determination being prevented by the uncertainties in (i), (iii) and 
(v) in $2. 

3.2. Nuclear Coupling 

The interaction between x-mesons and nucleons may be treated in a manner 
similar to that between leptons and nucleons except for a difference in coupling 
ea, P(N)=>N(P) +t) $H® cea (2) 
From (vi) §2, it appears that the interaction is quite strong, and may thus be an 
important source of nuclear forces at short distances. At high energies, therefore, 
the «-meson behaves in a manner similar to the 7-meson: the cross section for 
their production in nucleon-nucleon collisions is large, and the probability of the 
production of secondary x-mesons by collisions between incident «-mesons and 
nucleons is considerable. ; 

‘The strong interaction for the process (2) results in the production of many 
«-mesons in the cosmic radiation in the upper atmosphere. ‘The chargedx-mesons 
decay quickly into -mesons via process (1), and this process competes with the 
usually considered origin of the u-mesons in the decay of 7-mesons. From a 
study of the underground radiation, Bollinger (1951) has shown the latter source to 
be predominant up to energies of the order of 10" ev. This suggests, therefore, 
that the production of «-mesons in nucleon-nucleon collisions may be ten or more 
times less than that of 7-mesons at high energies. | Furthermore, the nuclear 
interaction of «-mesons will be smaller than that of 7—mesons by a similar factor. 
In what follows, an interaction length of a few tens of metres of water will be 
assumed for x-mesons. 


§4. POSSIBLE EFFECTS IN COSMIC RAYS 

Whereas most of the charged «-mesons produced in the atmosphere decay into 
#i-mesons on account of their short life, the neutral, «°-mesons produced either 
directly in nucleon-nucleon collisions or resulting from the decay of charged 
«-mesons can survive down to ground level on account of their lack of ionization, 
their fairly long interaction length, and also presumably of their lack of decay. 
The frequency of «°-mesons at sea level is then estimated to be as large as that of 
p.+-mesons. 

The «°-mesons which penetrate into the earth can produce charged «-mesons 


through the process «9+ P(N)>N(P)+«+(-), 22") 


if their energy is sufficiently high. ‘The charged x-mesons thus produced are 
lost mainly by decay and nuclear absorption, the ionization loss being negligible 
at the moderate depths and high energies in which we are interested. ‘The author 
has solved the diffusion equations for the charged «-mesons and their p-meson 
decay products under simplified assumptions. As is usual in such problems, a 
critical energy is defined, determined by the lifetime and interaction length of the 
«*-mesons, which is found to be of order 104ev. Below the critical energy, 
decay is predominant, and the x*-mesons are in equilibrium with the parent 
«°-mesons, and their intensity is small. On the other hand they supply a large 
intensity of ~*-mesons which increases up to a certain depth and then reaches 
saturation. Above the critical energy, the intensity of «*-mesons follows a 
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characteristic transition curve, possessing a maximum at a depth equal to the 
absorption length. At the maximum, the intensity of «*-mesons is found to be 
as high as that of the .*-mesons of the same energy. All effects resulting from 
the x°-mesons striking the ground are found to be negligible at depths greater 
than say 200 m water. ‘ 

The above picture suggests that there may well be an additional source of 
penetrating particles at moderate depths below ground, consisting of the charged 
x-mesons and their p-meson decay products. If this is so, at intermediate depths 
the depth-intensity curve should lie above that calculated from the energy 
spectrum of the 4.+-mesons at sea level, coinciding with it at depths greater than 
about 200m of water. Indeed, the experimental curve shows a gradual change 
of slope (Rossi 1948) while the momentum spectrum of the .~-mesons at sea level 
is well represented by a simple power law expression between 3 x 10° and 
3x10! ev/c. Although the measurements have not great accuracy, the 
difference between the two curves seems to be significant, and is very difficult 
to understand on the basis of our commonly accepted ideas. 

The appreciable number of high-energy «*-mesons expected to exist should 
give rise to nuclear interactions underground. ‘The first order process is similar 
to (2’) and would be observed as the annihilation of «*-mesons in flight. The 
second order process 
tc (NP) = NOP) Re) eR Oe a eee (3 
would be observed as the nuclear scatter of a fast particle and the third order 
process 


K+(k-) + N(P)—> P(N) +o +4-(kt) +O tnt) +H, owe. (4) 


would be observed as the production of an additional penetrating particle.* 
The cross sections for these processes become smaller by a factor of about ten 
at each step. 

These nuclear interactions may be contrasted to the nuclear interactions of 
p=-mesons, which latter are thought to be responsible for most of the nuclear 
disintegrations observed underground by George and Evans (1950) and Cocconi 
and Cocconi-Tongiorgi (1951). In the case of these p-meson induced dis- 
integrations, the mechanism appears to be first the production of a 7-meson in 
an encounter between the incident u-meson and a nuclear nucleon, followed 
by the re-absorption of the 7-meson inside the same nucleus (Hayakawa 1951). 
In the case of the «-meson interactions, the nuclear interactions of the secondary 
k-mesons are not strong enough to produce further nuclear collisions in the 
nucleus in which they are produced. Furthermore, on account of this weak 
interaction, the nucleon recoil will not be large and the processes under 
consideration should occur mainly without associated nuclear disintegrations. 

The processes (3) and (4) are suggested as a possible interpretation of the 
underground observations of scattering and penetrating pairs mentioned in the 
introduction. ‘The cross sections resulting from a detailed calculation are 
consistent with those observed, as also are the fairly frequent cases of consecutive 
scattering in each of two lead plates, and the small angular divergence of the pairs. 
Without some such speculation these phenomena are very difficult to explain. 


* I am indebted to Professor Rosenfeld for pointing out the possibility of the following alternative 
scheme for the production of a penetrating secondary particle : 


K*(K~)-+- N(P)P(N) + 9+ eo (et) + ot ata), 
The cross section for this process is found to be smaller than that of scheme (4). 
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The above picture is admittedly tentative, and only points favourable to the 
model have been discussed: it is quite easy to argue against it just as forcibly. 
The author’s intention has been to stimulate further experiments which might 
confirm or disprove the speculation. The experiments are not difficult and 
some are repetitions of old ones in which greater accuracy and reliability is 
needed. Possible experiments are as follows: 

(i) From reaction (2’), the frequent annihilation of charged particles in flight, 
and the frequent production of charged by neutral particles should be 
observed. This could be checked by a simple anticoincidence arrangement. 

(iu) An accurate re-determination of the intensity-depth curve, preferably in 
water, in order to decide whether this curve shows significant departures 
from that inferred from the measured energy spectrum of penetrating 
particles at sea level would be useful. 

(iii) According to processes (3) and (4), the particles involved in anomalous scatter 
and in penetrating pair production should all be heavy mesons, of mass 
of the order of 1000m,. Observations with a cloud chamber in a magnetic 
field might be useful in this respect. 
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Experiments on Production of Auroral Radiation 
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ABSTRACT. Further experiments have been made on the production in the laboratory 
of the ‘ forbidden’ red and green auroral lines, and the ultra-violet transauroral line of O1. 
The mean concentration of oxygen atoms in the 11S, and 5°D states in the columns has been 
found from absolute intensity measurements. The decay of the green radiation produced 
by a Tesla discharge through a high-pressure source has been investigated with a 
photo-multiplier tube. It occurs at approximately the rate to be expected if the atoms 
in the 1S state are undergoing spontaneous transitions to lower atomic levels, but 
complicating factors make it difficult to deduce transition probabilities from the experiment. 
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§1. INTRODUCTION 


HE principal ‘forbidden’ lines resulting from transitions between low 

| terms of O1, most of which appear in the spectrum of the night sky and 
aurora, have been repeatedly produced in the laboratory. Much still 
remains obscure, however, about the conditions under which they are excited. 
This paper contains an account of some quantitative experiments which have been 
done with discharge tubes containing gas at low pressure (§§ 2, 3) and atmospheric 


Forbidden ‘Transitions in O1 


‘Transition X (A) A (sec™) 
3P;—15o 2972°3 0-090 
116, 557768 2-20 
5P,—1D, 6300-2 0:0078 
5P,—!D, 6363°9 0-0026 


pressure (§4) with the aim of determining the excitation conditions more 
precisely. ‘The approximate wavelengths of the observed lines, and the transition 
probabilities (A) calculated from the electron configurations (Pasternack 1940), 
are shown in the table. 


§2. DISCHARGE TUBES 


The evidence on which design of low-pressure discharge tubes for production 
of the forbidden transitions can be based appears somewhat contradictory. Our 
earlier investigations (Emeléus, Sloane and Cathcart 1939) and those of Kvifte 
and Vegard (1947) showed that strong excitation was likely to occur in positive 
column tubes, if they were wide and the gas pressure relatively high, to hinder loss 
of metastable atoms by diffusion to the wall. Herman and Herman (1944) 
have, however, obtained both the green and red lines readily in tubes of only 
lcm radius passing feeble currents, whilst Hopfield (1931) has obtained them 
with narrow tubes passing heavy currents. Probe analyses which might account 
for these differences are lacking, and likely to be difficult to obtain because of the 
occurrence of irregular ionic oscillations (Sloane and Minnis 1937). The tubes 
made for the present investigation have been larger versions of the 7-shaped 
tubes used previously (Emeléus, Sloane and Cathcart 1939), with large hollow 
iron electrodes in the side arms and horizontal sections up to 6cm in diameter 
and 110cm long. A graded quartz seal was recessed longitudinally into one end 
of the horizontal section for observation of the ultra-violet line. In our earlier 
investigations it was found that the green auroral line could only be obtained when 
the positive column had irregular moving striations. As there appeared, therefore, 
to be little advantage in use of d.c., we have used 50 c/s a.c. at between 1500 and 
4000v. Moving striations were present in each half cycle. Most experiments 
were done with the oxygen diluted with argon, and a few with oxygen and neon. 
The former gave the better source. ‘The forbidden lines were not emitted from 
the negative glow and Faraday dark space. 

The tubes behaved in most respects as had been anticipated and the green 
line was obtained consistently. Using argon, it was also found, in agreement 
with our earlier results, that it appeared almost in isolation at the periphery of 
the discharge and where the discharge diffused into the corner at the junction of 
the horizontal and vertical sections of the tubes. The dipole lines of Or and 
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AI were much more concentrated near the axis of the discharge. Similar, but 
less marked, radial distributions were found with oxygen and neon. An un- 
expected advantage of use of a.c. was that there was less tendency for the oxygen 
to move away from the middle section of the length of the tube than with d.c. ; 
this may be an electrophoretic effect. It was again noticed that the line could 
sometimes be obtained more brilliantly after the tube had been open to the air, 
and only torched before admitting gas, than when it had been degassed rigorously, 
possibly because imperfect degassing leaves surface films on the glass which 
hinder recombination of oxygen atoms. The red lines were also obtained, but 
less consistently, probably because the long life of the 1D state (see table) makes it 
necessary to operate within a restricted range of conditions in which the meta- 
stable atoms are not too rapidly destroyed by collisions with electrons or molecules ; 
they were not observed preferentially from the periphery of the discharge in the 
same way as the green line. ‘The transauroral line was obtained without difficulty, 
although not with high intensity; the measured wavelength was 2972:32A 
(Sayers and Emeleus 1950). When the auroral lines were excited in argon— 
oxygen discharges, a patch of continuous spectrum usually appeared between 
about 6400 and 6550A. It was always present when the red lines were obtained. 
Its origin is unknown. 


§3. INTENSITY OF THE GREEN AURORAL ‘LINE 
AND OTHER SPECTRA 


If emission of a line occurs spontaneously, the rate EH of radiation of energy 
of that wavelength per cm? of the discharge tube will be NAhverg/sec, where 
N is the concentration of excited atoms, v the frequency and A the transition 
probability. If EF can be measured and A is known, N can be calculated, 
neglecting self-absorption in the discharge, which is usually permissible. ‘This 
provides the most direct method at present available for finding concentrations 
of excited atoms. It has been used for part of the green spectrum from argon— 
oxygen discharges by comparing the intensity of the radiation received at the 
slit of a Hilger constant deviation spectroscope with that received from a calibrated 
hot-filament lamp. One half of the length of the slit was illuminated by light from 
the middle of the discharge tube, viewed laterally, and the other half by light 
from the standard lamp, both being brought to the slit by nearly identical optical 
systems. ‘The intensity of the green light from the two sources was generally 
matched visually, but in a few cases the spectra were photographed and micro- 
photometered. It was found convenient for visual observation to use a slit 
0-3 mm wide. The lamp then gave a brighter illumination than the discharge, 
and the light from the former was therefore reduced by ground-glass plates 
calibrated for absorption. Since the discharge was intermittent, and inhomo- 
geneous radially, more accurate photometric measurements were not warranted. 

Some results are given in figs. 1(a) and 1(6). Before making each measure- 
ment a fresh charge of gas was admitted. Fig. 1(a) shows on a logarithmic 
scale the number of quanta J emitted per cm?/sec, as a function of tube current 7, 
for a tube 4cm in diameter containing 0-4mm oxygen and 3:0mm argon. ‘The 
lines 1, 2, 3 in the figure are for the auroral line (5577 A), the unresolved permitted 
Or multiplet at 53294 (35P —5°D, excitation potential 13-0v), and for a region 
of 204 centred on 55304 of the continuous background which appears from this 
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source (Emeléus, Sloane and Cathcart 1939). Allowance has been made for 
the continuum in calculating the intensities of the lines. The three curves 
are nearly parallel. For the larger currents [x 7i*. Kvifte and Vegard (1947), 
using mixtures ot oxygen and neon, found that [7% for the auroral line and 
the permitted line at 5555A (3°P —7°S, excitation potential 13-2 v). 

Fig. 1(b) shows J for the same spectra as a function of partial pressure p of 
oxygen, in mmHg, for the same tube passing an r.m.s. current of 0-10 amp at 
a total pressure of 4:0mm of oxygen and argon. For p<1-7mm Jax p*’, 
and for p>1:7mm I «<p-*1. The change in slope is connected with a change 
in type of the positive column; over the whole pressure range there was marked 
radial constriction of the discharge, but above 1-7mm this increased more 
rapidly with increase of pressure than below 1:7 mm. 
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Fig. 1. (a) Mean rate of emission of quanta J plotted against current 7. (b) Mean rate of emission 
of quanta J plotted against pressure of oxygen p. 


The greatest value of N obtained for the 1S state was 1-5 x 10"/cm’. The 
corresponding concentration of oxygen atoms in the discharge, assuming all 
the oxygen was dissociated, was 4 x 10!”/cm’. ‘This is so much in excess of the 
number of excited atoms that the results, even allowing for very incomplete 
dissociation, are consistent with spontaneous emission of the auroral line. A 
previous estimate made by comparison of the intensities of the auroral line and 
an AI line gave a concentration of excited atoms of 2:5 x 10"/cm® under similar 
conditions (Emeléus, Sloane and Cathcart 1939). The concentrations of 
oxygen atoms in the 5°D states are of the order of 10-7 or 10-8 times smaller, as the 
intensities of the lines at 5577 A and 5329A are comparable, and the values of A 
for the permitted lines are of the order of 107-108 sec~!. No information can be 
deduced about the emitter of the continuum, whose origin is uncertain. 
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Analysis of the data shown in fig. 1, to find the probability of de-activation 
of excited oxygen atoms on collision with oxygen molecules and inert gas atoms, 
by a method developed by Kvifte and Vegard (1947) for neon, gave negative 
values for the probabilities. This might be due to regeneration of excited 
oxygen atoms occurring in argon by some process not operative in neon. We 
think it is more probably due to the constriction of the discharge in argon, which 
we found to be greater than in neon, making invalid Kvifte and Vegard’s 
assumption that the current density is uniform over the cross section, and the 
consequences of this. A quantitative diffusion theory similar to Kvifte and 
Vegard’s could be developed for the region between the core of the positive 
column and the wall, with a solution differing mainly from Kvifte and Vegard’s 
in introducing Bessel functions of both kinds. There is, however, no point in 
doing this, as the boundary conditions at the edge of the core, which is itself 
diffuse, are uncertain. Caution is also necessary in using Kvifte and Vegard’s 
results for the de-activation coefficients of the 1S state of O in O,—Ne mixtures, 
as similar results could be obtained if de-activating collisions occurred with high 
efficiency, with impurities, or ozone, present in minute quantities. 


§4. HIGH PRESSURE SOURCES 


Although the discharge changes from a glow to a spark as the pressure is 
increased from a low value to atmospheric, it has been shown by Janin (1946) 
that it may be a good source of forbidden spectra. Vegard and Kvifte (1948) 
and Jenkins, Bowtell and Strong (1949) have found that the afterglow of Tesla 
discharges at or near atmospheric pressure through a mixture of argon, or argon 
and nitrogen, with a trace of oxygen of the order of one part in 10°, is a good source 
of the green auroral line, which is accompanied by diffuse bands and a continuum. 
This gas-filling is used for certain commercial incandescent lamps. We have 
made some studies of this interesting source, using 300 and 500 watt Osram bulbs 
and the tube shown in fig. 2 with nickel rod electrodes. ‘The general features of 


Fig. 2. High pressure discharge tube. Length 21 cm. 


the discharges were the same in all cases. Near the point of entry to the gas, 
from the metal electrode or through the glass, there is a local spark, which 
branches at greater distances, and is surrounded by a green aureole a few 
millimetres in diameter. Green clouds move upwards from the spark, presumably 
on account of thermal convection of excited gas. ‘These spread and, if the dis- 
charge is stopped, mingle with a general green afterglow which has been 
detected visually for more than 5 seconds. ‘I'he spectrum of the spark contains 
nitrogen bands; these are weak or absent from the spectra of the aureole around 
the spark and from such afterglows as we have examined. Cyanogen bands also 
appear erratically, usually with feeble intensity. ‘The main feature of the 
spectrum of the aureoles and the afterglows is the green auroral line. The 
continuum and bands near 5577 A were most marked in the spectra of the aurecles, 
and less strong, if present at all, in the spectra of the sparks. A microphotometer 
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trace of the spectrum of an aureole is shown in fig. 3. The peaks of intensity 
at 5584, 5570, 5550, 5530, 5490 and 54554 can be attributed to an unstable 
AO molecule (Herman, Weniger and Herman 1951). 

The decay of the green light from the afterglow has been followed by a photo- 
multiplier, the current from which was amplified and shown on a cathode-ray 
oscillograph. For about the first 0-2sec after the exciting discharge was 
discontinued the decay of light from any selected part of the tube was irregular, 
in some instances actually increasing with time for short intervals. After about 


eH head 


5450 5500 i 5550 5600 
Fig. 3. Microphotometer trace of auroral line and associated band for high pressure source. 


0-2sec until about 1-5sec, after which the light was too feeble to be recorded 

in this way, the decay was exponential. The values of the decay constant were 

in the neighbourhood of 2 to 3 sec"!, depending on the tube used and the part 

studied. The following data show the trend of the values of the decay constant 

and the consistency obtained. 

(i) First 300 watt Osram lamp: discharge scanned in middle of bulb; decay 

constants; 2°5, 2°6, 2:3, 2:9, 2:9 sec}. 

(ii) Second 300 watt Osram lamp: discharge scanned in middle of bulb; decay 
constants, 2:4, 2:2, 2:2sec!. 

(111) Electrode tube (fig. 2) : discharge scanned close to electrodes ; decay constants, 
Soe SCE. 

(iv) Electrode tube (fig. 2): discharge scanned at end remote from electrodes; 
decay constants,~2-6, 2:4, 2°7, 2'6; 2:9 sec *. 

The decay constant was usually greatest close to the exciting discharge. 
Further away it usually diminished, but was not found to fall below about 2:0 sect 
in these experiments. 

At least four processes may be responsible for the decay: (a) convection of 
the luminous gas; (b) diffusion to the wall of excited atoms in the 1S state and of 
the associated excited AO molecules; (c) emission of radiation from the same 
atoms and molecules; (d) de-activation of the excited atoms and molecules in 
collisions of the second kind. It is not possible to distinguish certainly between 
them. It has been shown that (b) can lead to about the observed rate after the 
effects of initial irregularities due to localization of the discharge and convection 
have died down (Emeléus, Sayers and Bailey 1950). On the other hand, the 
least value found for the decay constant is close to the A coefficient for the green 
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line (see table). More recently, however, L.. Herman (private communication) 
has found values of the decay constant of 1-5sec-!. If the value of A given in 
the table is accurate, this indicates, as he points out, that regeneration of oxygen 
atoms in the 'S state may be taking place. 


§5. DISCUSSION 


The experimental problem of finding sources of the forbidden spectra 
coming from transitions between low terms of O1 (and similar spectra of other 
light neutral atoms and molecules) is now largely solved. New methods for 
producing the necessary concentrations of metastable atoms have been developed, 
using gas at high pressure, which supplement the original methods employing 
gas at low pressure, and the original low pressure sources have been further 
developed. However, there is not yet available a satisfactory detailed description 
of the condition of the gas in any of the sources, for the criteria for a discharge, 
to serve as a good spectroscopic source for this purpose and to be amenable to 
electrical analysis, remain incompatible. The experiments made since our earlier 
analysis (Emeléus, Sloane and Cathcart 1939) have apparently all been carried 
out with discharges that were energized by an alternating or intermittent voltage. 
In some experiments, such as those at low pressure described in the present 
paper, a.c. has been used deliberately. With gas at high pressure it is probable 
that the discharge is inherently intermittent, wandering from one part of the gas 
to another if it is not interrupted in a single locality. ‘The new experiments at 
low pressures, and the experiments which have been done at high pressure, do 
nevertheless confirm our earlier conclusions that forbidden transitions are likely 
to be produced in most isolation in feeble discharges and in’the aureoles round 
heavy discharges or in afterglows, i.e. under conditions generally similar to those 
in the upper air, and that the metastable atoms may be produced by more than 
one process. 

It is still uncertain if forbidden lines can be produced with appreciable 
intensity from the plasma of a heavy current discharge. Relations such as those 
shown in fig. 1 cannot be interpreted definitely without more detailed knowledge 
of the condition of the ionized gas than is yet available. It is not unusual for 
simple empirical formulae to be found to hold approximately for discharges under 
circumstances which indicate that the apparent simplicity is due to mutual 
compensation of a number of complex conflicting processes. ‘The irregular 
moving striations present in the tubes with d.c. supply, which automatically 
introduce afterglows, might be suppressed by building an anode with a negative 
anode drop in potential (Pupp 1933), if this would operate steadily with oxygen 
present. Irregular low-frequency fluctuations may however occur even in the 
absence of moving striations (Cobine and Gallagher 1947), so that elimination 
of moving striations would not necessarily leave the discharge steady. ‘The 
existence of striations and fluctuations, and even the use of 50 c/s a.c., does not 
nevertheless completely invalidate the application of the quantitative theory 
of de-activation developed by Kvifte and Vegard (1947), since the lifetime of 
the metastable atoms is much greater than the reciprocal frequency of the greater 
part of the electrical noise spectrum. ‘The comparison of de-activation times 
with lifetimes of metastable atoms is simpler in afterglows where applied electric 
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fields are absent than when a current is passing, but even with afterglows re- 
generation of metastable atoms may add to the complexity of the processes (§ 4), 
whilst, as with sustained discharges, the determination of the degree of molecular 
dissociation and of the concentration of free radicles presents difficult problems 
not yet satisfactorily solved. 
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LETTERS: TO PHE* EDITOR 


Helium II Transfer on Metal Surfaces 


The experiments of Mendelssohn and White (1950 a, b) on the transfer rate of the 
He II film had shown that, whereas clean glass surfaces yielded values almost identical 
with the earlier work (Daunt and Mendelssohn 1939), metal surfaces always showed much 
higher transfer rates. "These results were confirmed by the observations of Dash and 
Boorse (1951). ‘Two explanations of this phenomenon seemed feasible, either that the 
different transfer rates may be due to the chemical nature of the substrate, or that there 
may be only one transfer rate (that exhibited by the glass surfaces) and any higher rates 
would be due to the micro-structure of the substrate. For a number of reasons we had 
preferred the second explanation, and in particular considered that narrow surface cracks 
of the order of magnitude of the film thickness or somewhat larger might be filled up with 
liquid helium and thereby produce a much greater transport of liquid. 

An opportunity to test this hypothesis arose when the Metrology Division of the 
National Physical Laboratory offered to prepare for us stainless steel beakers of a much 
better finish than any metal surface that had been previously used for helium transfer. 
The transfer rate obtained on two beakers of this kind in dependence on temperature is 
shown in the accompanying figure. As can be seen, the results for both beakers are 
identical, and the transfer rate is practically the same as that for a baked-out glass surface. 
Subsequently, one of the beakers was brought to red heat, thus destroying the smooth 
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surface finish, and it was then lightly polished with cotton wool to remove any deposit 
of oxide. Examination under a low power microscope did now show cracks and draw-marks. 
A transfer experiment on this beaker showed indeed that the flow rate had changed 
considerably, being now more than three times the original value. 
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Transfer rates of the helium II film on two highly-polished stainless steel beakers (O and @ ), and 
on one of these beakers (©) after destruction of the surface finish by heating (@). The 
broken curve is that obtained for baked-out glass by Mendelssohn and White. 


In conjunction with other flow experiments and observations on the formation of the 
helium film which will be reported in detail later, we conclude that the present experiments 
show that the high transfer rates observed are in fact due to the micro-structure of the 
substrate. 

We wish to thank Messrs. F. H. Rolt, V. W. Stanley and S. Woodhouse of the Metrology 

Division, National Physical Laboratory, for their helpful co-operation. 


The Clarendon Laboratory, B. S. CHANDRASEKHAR. 
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Anomalous Thermal Conductivity of Pure Metals at 
Low Temperatures 


In a recent paper on the thermal and electrical conductivities of sodium at low 
temperatures, Berman and Macdonald (1951) have drawn attention to two types of 
discrepancy between their experimental data and the behaviour predicted by the electron 
theory of metals. They observe firstly that the thermal conductivity fails to pass through 
a theoretically predicted minimum at a temperature in the neighbourhood of 0:25 ©, where 
© is the Debye temperature, and secondly that, although the experimental Lorentz number 
K/oT (where K and o are the thermal and electrical conductivities respectively) exhibits 
a minimum at low temperatures, in accord with theory, the location and shape of the 
minimum deviate considerably from the theoretical results of Makinson (1938). The 
purpose of this letter is to point out that the discrepancy for sodium at the lowest 
temperatures (< 0-1) is essentially of the same type as that already well known for ten 
other metals studied by Griineisen and Goens (1927), de Haas and Rademakers (1940), 
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Hulm (1950) and Webber, Andrews and Spohr (1951), and originates in the fact that the 
observed coefficient of the T? term in the electronic thermal resistivity is nearly an order 
of magnitude lower than the theoretical value. A placing of special emphasis on the Lorenz 
number for sodium tends to obscure the fact that it is in the thermal resistivity rather than 
in the electrical resistivity that the main discrepancy occurs. 

Makinson’s theory leads to the result that for temperatures below 0:10 the electronic 
thermal resistivity may be expressed as the sum of a lattice vibration scattering term aT” 
and an impurity scattering term 8/7. The coefficients may be put into the form 


_ 95:3N,2" 
7 Ree 


where N,, is the effective number of conduction electrons per atom, K., is the electronic 
thermal conductivity at high temperatures, py is the residual electrical resistivity and Ly is the 
Lorenz constant }(7k/e)?. Assuming one electron per atom and substituting the other 
appropriate constants for sodium one obtains an « value of 2:26 x 10~* watt/cm/deg*, which 
is about six times Berman and MacDonald’s experimental figure of 3:8 x 10~* watt/cm/deg’. 
Alternatively, substitution of the experimental a into eqns. (1) leads to a value for N, of 
0-069, which is of the same order of magnitude as the values for N, obtained for Sn, Hg, 
Pb, In, Cu, Au, Al, W, Fe and Pt by the same method (see Hulm 1950, table 2), and is 
undoubtedly much too low. It should be mentioned that recent refinements of Makinson’s 
theory proposed by Umeda and his co-workers (1941, 1949) and Sondheimer (1950) have 
decreased the theoretical value of « by about 25°, which is, however, far from sufficient 
to remove the disagreement with experiment. 

While anomalous heat conduction in metals at temperatures above 0-10 may in part 
be due to deviations of the phonon spectrum from a Debye pattern, as Berman and 
MacDonald have stated, the quantitative discrepancy in the lattice vibration scattering 
term of the electronic thermal resistivity below 0:1 evidently involves a direct breakdown 
of the free-electron theory. The situation that an ‘ ideal’ metal such as sodium gives a 
departure from theory which is of the same order of magnitude as those for more complex 
metals such as Sn, Pb, W, Fe and Pt suggests that this breakdown is of a rather fundamental 
nature. J. IK. SEhoune 


Institute for the Study of Metals, 
The University of Chicago. 
12th December 1951. 
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Transition Probabilities of some Band Systems of Nitrogen 


With a view to understanding the processes in the upper atmospheric sources of the 
First Positive, Second Positive and Negative band systems of nitrogen, Bates (1949) 
derived the relative probabilities of transitions between the various vibrational levels by 


computing the square of the overlap integral containing the normalized vibrational wave 
functions, 


(7 2 
| Rin’, v'/7) Rn", v" |r)? dr\. 


Using Morse’s potential-energy expression, Bates solved the relevant differential equation 
in terms of Laguerre polynomials. 
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In verifying his theoretical predictions for the Second Positive system, Bates used 
the extensive data derived by Tawde (1934) from laboratory measurements of peak 
intensities. He pointed out the limitations of these data (due to wide differences between 
values from one condition of excitation to the other) and the importance of accurate 
quantitative intensity measurements. It would appear, however, that he was unaware 
of 'Tawde and Patankar’s (1943) precise experimental intensity observations, the first to 
be computed from the areas under the contours of intensity curves. These observations 
should give a better approach to reality since they account for any variations in the overall 
intensities in a band. They are used here to test Bates’ derivations, which, at the same 
time, are compared with the transition probabilities already computed (by Langstroth 
(1935) and by Tawde and Patankar, independently) from Hutchisson’s (1930) intensity 
integral and with those calculated by Nicholls (1950) using the method of Gaydon and 
Pearse (1939) and Pillow (1949). The results for the Second Positive system are collected 
in the following table, from which it will be seen that Tawde and Patankar’s revised observa- 


tions show much better overall agreement with Bates’ derivations than do Tawde’s earlier 
data. 


Band Experimental values Theoretical values 
ooo oF : ee SP ee a RI. 
Oe (1) (2) (3) (4) (5) (6) 
0,0 0-40 0-43 0-48 0-53 0-48 0-46 
Ord 0-36 0:32 0:32 0:27 0-38 0-38, 
Ox: 0:17 0-15, 0-14 0:12 0-12 0-12, 
O53 0-05, 0-07 0:04, 0:05 0-01, 0-02 
te) 0-49 0-44, 0:39 0-43 0:32 0:32, 
sla 0-02, 0:02, 0-03, 0-04 0-03; 0-02, 
Ae 0-16 0-16, 0-22 0-21 0-34 0-32, 
3 0-20 0-20, 0-19 0-12 02287 0-24; 
2,0 0:29 0-22, 0-12 0-21 0-15. O21 
Deol 0-33 0-49, 0:37 0:39 0°25, 0:30 
Ded. 0-03 0-03, 0:00, 0-01 0-02 0-04 
Das 0-05, 0-06, 0:09, 0-10 0:22, 0-25 
2,4 0-14 0-17, — 0-08 0:34 0-32 


(1) T'awde (mean values); (2) 'Tawde and Patankar; (3) Bates; (4) Nicholls; (5) Tawde and 
Patankar (Hutchisson’s method); (6) Langstroth (Hutchisson’s method). 


As a consequence of a number, of simplifying assumptions, as well as the different 
values of molecular parameters used in the calculations, there is a variation amongst the 
different sets of theoretical values. The sensitiveness of values of transition probabilities 
to a slight change in molecular constants has been demonstrated by McKellar and 'Tawde 
(1951). Proper comparison would require theoretical calculations from the same molecular 
constants; with the limitations of the available data, however, the following broad features 
become apparent. 

Bates’ and Nicholls’ values do not differ widely, except for the (1,3) and (2,0) bands. 

The two sets of values derived from Hutchisson’s theory are in a class by themselves, 
comparable with one another and departing about equally from Bates’ or Nicholls’ results. 
This is to be expected since they follow from the same theoretical expression, true only 
' for a harmonic oscillator; the small differences between the two arise from the use of 
different molecular constants in the two sets of calculations. All four sets of theoretical 
results agree fairly well amongst themselves for bands associated with lower quanta. 

In view of the errors and uncertainties in photographic photometry and the limitations 
of theoretical results, there is remarkable agreement of experimental observations (particularly 
those of Tawde and Patankar) for many bands, even bands of higher quanta, with Bates’ 
or Nicholls’ values. This relatively closer agreement over the entire region is attributable 
to the fact that, in different ways, both take account of anharmonicity: Bates directly 
uses the experimentally justified Morse potential-energy expression and normalized wave 
functions in the integration, while Nicholls starts with the parabolic curve of a harmonic 
oscillator and distorts it to fit the Morse expression. As Hutchisson’s theory is based on 
a harmonic oscillator only, it is at higher quanta, e.g. in the (0, 3), (2,3) and (2, 4) 
bands, that the values in the last two columns of the table show greater disagreement 
with Bates’, Nicholls’ and the experimental values. 
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More work on these aspects of theories of transition probabilities in relation to 
experiment is in progress. 
We are grateful to Dr. N. R. Tawde for suggesting this study and for very helpful 
discussion in the preparation of this note and to Dr. W. Jevons in England for kindly 
revising the manuscript to the size of a ‘‘ Letter”’. 
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ABSTRACTS FOR SECTION B 


Characteristics of Radio Echoes from Meteor Trails: III—The Behaviour of the 
Electron Trails after Formation, by J. S. GREENHOW. 


ABSTRACT. This paper deals with the durations and amplitudes of the radio echoes 
reflected from meteor trails. For wavelengths of 4:2 and 8:-4m the durations of most 
echoes are shown to be very nearly proportional to the square of the wavelength, as predicted 
by simple diffusion theory. It is suggested that long duration meteor echoes are due to 
total reflection from trails of very high electron density, and it is shown that the amplitudes of 
these echoes would only be of the same order as the amplitudes of the short duration echoes 
from trails of low electron density discussed in an earlier paper. Further evidence is given 
to show that the amplitude fluctuations.in the radio echoes are caused by distortion of the 
column of ionization, relative motion of several reflecting centres due to wind turbulence 
producing changing interference conditions. Turbulent winds of the order of 20 m/sec are 
inferred at altitudes of 80 to 100 km. 


The Propagation of Sound through Gases contained in Narrow Tubes, by 
L. E. LAWLEY. 


ABSTRACT. A comprehensive investigation is made into the validity of the Helmholtz— 
Kirchhoff equations for the absorption and velocity of sound in a gas contained in a narrow 
tube. Results are given for air, oxygen, nitrogen and hydrogen contained in tubes of radii 
between 0:15 mm and 1:17 mm at frequencies between 60 and 150 kc/s. The form of both 
the velocity and absorption equations is confirmed, but the viscosity-thermal conductivity 
constant is found experimentally to be about 5° above the theoretical value. A method is 
outlined in which the free space absorption can be calculated from measurements in a tube 
at varying pressures, and results are given in the case of oxygen. 


The Magnetic Electron Microscope Objective Lens of Lowest Chromatic Aberration, 
by G. LIEBMANN. 


ABSTRACT. ‘The shortest focal length f and smallest chromatic aberration constant C, 
of an objective lens is directly proportional to the square root of the (corrected) accelerating 
voltage V, and inversely proportional to the maximum magnetic field strength H,, in the 
parallel part of the pole-piece gap. For given values of V, and H,, both f and C, fall 
monotonically with the relative pole-piece separation S/D, although the reduction is insig- 
nificant for S/D >2. For each lens geometry there exists an optimum set of lens dimensions. 


Effect of Strain Ageing upon the Distribution of Glide Bands in Cadmium Ces, 
by D. F. GrBBons. 


ABSTRACT. It is shown that there is a marked difference in the distribution of 
visible glide bands, after plastic deformation, bevween cadmium crystals containing 
small amounts ot nitrogen and those relatively free from nitrogen. The distribution of 
glide bands in the crystals containing nitrogen is compared with Liiders bands in strain- 

aged mild steel. ; 


Single Contact Lead Telluride Photocells, by A. F. Gipson. 


ABSTRACT. In a recent paper the photoconductivity of lead sulphide films was 
discussed in terms of a barrier mechanism. Large single crystals of lead telluride have now 
been grown and, with these, single barrier photocells can be prepared. The simplest barrier 
is formed at a metal-semiconductor contact. Such barriers, as exhibited in rectifiers, have 
been mainly studied. It is found that the photosensitivity and time constant of such cells is 
determined by the total current at the contact point, as expected from barrier theory. The 
yariation in barrier height with temperature and time has also been determined, 
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Lead Sulphide Rectifier Photocells, by A. F. GIBSON. 


ABSTRACT. Following the study of lead telluride rectifier photocells described in a 
contemporary paper, similar properties have been found for lead sulphide. 


A Note on the Structure of Selenium, by E. BILLic. 


ABSTRACT. The recently established identity of the width of the (optical) energy 
gap, 2-05 ev, for both amorphous and crystalline selenium, suggests that the chain 
structure of crystalline selenium is maintained also in the amorphous phase. The low 
electrical conductivity of amorphous selenium is ascribed to potential barriers existing 
between individual chains. 


Electronic Structures and Physical Properties in the Alloy Systems Nickel-Copper 
and Palladium—Silver, by B. R. COLEs. 


ABSTRACT. A comparison is made of the magnetic, electrical and thermo-electric 
properties of alloys in the two systems, and the lattice spacings, optical properties, electronic 
specific heat, and x-ray spectroscopy of the nickel—copper alloys are also considered. 
Palladium-—silver alloys show the behaviour expected of them on the band theory of metals, 
but nickel—copper alloys show well-marked and consistent deviations. It is concluded that 
unoccupied d-electron levels exist in the copper-rich alloys, whereas such levels do not 
exist in the silver-rich palladium-silver alloys. From theoretical considerations it seems 
unlikely that the unoccupied levels form part of a collective d-band, and alternative models 
are discussed. That the two systems differ so markedly is surprising, but a tentative 
explanation of the differences is sought in the differing proportions of the atomic volume 
occupied by the ions in copper and in silver. 


Torque Curves and other Properties of Permanent Magnet Alloys, by K. HosELitz 
and M. McCaic. 


ABSTRACT. Further investigations of Alcomax III and related alloys have been made 
using a torque magnetometer. Samples with both columnar and random crystals have been 
used. ‘Torque curves have been obtained after various sequences of heat treatment designed 
to throw light on the changes that occur in these alloys, and the magnitude and direction of 
the magnetic field applied during cooling has been varied. A method of deriving three 
independent constants of crystal anisotropy without using single crystals has been devised. 


